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The Laplace cascade laethod is concerned with second order linear 
hyperbolic equations of the fora 



(1) + a(x,y) u^ 

The substitutions 
equations 



■b(x,y) 

» Uy + au 



u + c(x,y) u 0 . 
and u_^ * Ujj; 



lead to the 



( 2 ) 

and 



XL^ + bu^ - hu * 0 



(3) 



u, + auj^ - ku • 0 , 



where h * a + ab - c and k * b + ab - c are the two Darboux 
i y 

invariants. If either h or k vanishes, (1) has been reduced to a 
system of two first order equations, while if neither vanishes the 
equation may be cascaded in two directions. Solving (2) (or (3)) for 
u in terms of u^(or and siibstituting the resulting expression 

into (1) yields an equation for u^(u^jl) of the same form as (l) but 
with new coefficients, in general. This process may be Iterated, forming, 
a chain of equations, until either the original equation reappears, or one 
of the corresponding invariants vanishes. 

An extension of Volterra’s product Integral to the non-homogeneous 



system 



is aade first. Then tlie Laplace laetliod is exteiided to systems of 
second order hyperbolic equations, of the form 



^x^y J = 1 ^ X J*1 jy~ J=1 

, i = 1>2, ♦, «,n , 

Matrix notation is promptly Introduced, and equation (5) is revrltten 
as 



( 6 ) U + AU -t- BU -h CU - 0 , 

JL ^ 

The two related matrix invar iants for (6) are 

H « Ax + BA - C 
K = B 4- AB - C , 

and the chain of equations is developed as for the single equation. If 
either H or K is identically zero, the resulting two systems of first 
order equations, may he solved hy employing the ahove-mentloned extension of 
Yolterra’s product integral. 

The invariance of H and K are discussed in a manner analogous 
to the invariance of the functions h and k of the single equation. 

This is followed hy a consideratlan of periodic systems, l,e, systems siich 



r 




that after j Iterations the original equation reappears. This dis- 
cussion results in two theareas, the first of which is 



Theorea I - A system of equations of the fora (6) having 
constant matrix coefficients A and B, can he reduced to the form 



U_ 



t I 

H W 



Ey 



hy a change of variables U » , if an only if AB * BA. The second 
theorem, arising from systems of period two leads to a discussion of the 
form of the solution to the natrlx analog of Liouville’s equation 

’ ( 3 ) 




Discussion of the form of the solution when the chain terminates 
after a finite nuiiber of iterations leads to two further theorems which 
are completely analogous to theorems proved hy Darhoux for the single 
equation. 

The second ext^asion of the Laplace cascade method is made to the 
third order linear hyperbolic equation in three Independent variables of 
the form 



Here the number of invariant functions to be consider Juiiq)s from two to 
sixteen. The nature of the "invariance* of these functions is different 
from that of the second order invariants, in that some of the functions 
are time invariants, while others can only be considered as quasi -invariants. 
Four methods of cascading the equations are discussed, each of which requires 
severe restrictions on the coefficients a, b, ...,g. The class of equations. 



for which each method will result in a termination of the chain after a 
finite number of iterations, is explicitly pointed out. 

The final extension is a generalization of this third order exten- 
sion to the n^^ order linear hyperbolic equation with n independent 
variables, of the form 



( 8 ) 



.n 


















+ • • -tb-Ux + cu = O 



Introducing the linear operator D, so that (8) becomes 

D(u) = 0, 

the number of identities which may lead to a decomposition of (8) into 
a system of two equations of lesser order is computed. The exact number 
of true invariants is determined to be n(n-l), while upper and lower 
bounds on the number of quasi-invariants are forratilated. Theorem V proves 
that an invariant h is a true invariant , i .e . invariant under the change 
of coordinates u = ^(x 3 ^,X 2 >» u # if 



for some i,j = 1,2, ...,n, i j, vrfiile ar^ other invariant is, in general 
not a true invariant, hence a quasi-invariant. 

A discussion of earlier extensions is included. 
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SECTION I 



INTRODUCTION 



The cascade method was devised hy Pierre Simon Laplace^ and 
investigated in complete detail hy Gaston Darhoux.^ Earlier ex- 
tensions were made hy Uo Dini^3 and J, LeRoux^, and these exten- 



sions are discussed briefly in SECTION Y1 of this paper. The method 
itself deals with the linear hyperbolic equation with variable coef- 
ficients which has the form 



where the coefficients a^ b, and c are analytic in a certain domain 
Dj or at least as many times continuously differentiable as we feel is 
necessary, and u — u(x,y) is a real function of real variables. 

Laplace's cascade method begins with the Introduction of ^i^t 
are now called the two Darboux invariants, h and k, which are defined 
by the relations 



( 1 ) 




b(x,y) ^4- c(x,y)u = 0 



h — -f ab = c 



^x 



( 2 ) 



Ic =. ^b 

^y 



S. ab - c 



« 



Then (l) may be written in either of the following forms; 



1 






) 



c 



2 



( 3 ) 



^ + au) + b ( ^ + au) - hu = 0 j 

+ bu) + a( + bu) - ku = 0. 

07 ox Ox 



We may then consider the obvious substitutions 



(h) 




au , 



^ + bu , 



and further reduce (3) to 



(5) 



Jy 



+ bvi^ " hu = 0 j 



+ au - ku = 0 , 

"l 



Now if we should be so fortunate that either function h or k is 
identically zero, we will have succeeded in reducing our original 
second order equation to a system of two first order equations, 
either 



u 



-JIT 



1 = B u 

d 7 



+ au 



0 j 



(6a) 



1 

I 

‘u 



3 






or 



(6l>) 






bu 






< + au - =0 

^7 "1 -> 



either of which system may he solved hy quadratures. 

If, however, as is more likely the case neither h nor k 
is identically zero, all is not lost. To h© specific, let us consider 
the system (6a), as the details are quite similar whichever system 
we choose Instead of (6a), we have 



\ — + au j 

^7 ^ 



(6a«) 






-J- hu, - hu = 0 



We may integrate the first equation of (6a’) to obtain 



( 7 ) 



u = e 



© 



d7 X(x) 



where l(x) is an arbitrary fujaetion of x. Substitution of this ex- 
pression for u into the second equation of (6a’) yields 

jaaj 



^ U , 



4- 



he 



Cady 

dv +■ 



dy X(x) 



= 0 





% 




( 




h 



or 






ady 



( 8 ) 





+ bu 


. f 




1 


1 


, 


— 





e 



ady 



dy + X(x), 



Taking the partial derivative of both sides with respect to the vari- 
able y gives 
I ady 






L 


1 

^ + bu 


log h 


. bu' 




1 


a y 

‘ , 1 


ax ^ 



» v> t 3 b u 

by 3 



jady 



u 



After some straightforward algebraic manipulations > this becomes 

p 



^ X ^ y 



a ^ \ 
+ 1 1 



a ^ 



(9) 



1 

a y 



. c = 



0 



where a 



= a - 



a 



log 



^ 7 



b, and 



c- = c - 



^ a + 

3 ^ 



^ b 

■Jy 



^ log h 

'Ty 



We observe that (9) is of the very same form as (1), but with new 
coeffients a^, b^> and . Hence we may repeat our original pro- 
cess, in hopes that either of the nevx invariants h^ or will 

perhaps be zero. If not, we may proceed to evolve a new equation for 
the variable u^, and in fact we may iterate the entire procedure as 




■I 



i 



5 



often as necessary to produce a chain, or cascade of equations, in the 
hope that at soiae point the iteration will stop heeause one of the in- 
Tariants will heeoiBe zero. This will enahle us to solve a first order 
system hy quadratures, and hy tracing o\ir way hack through the chain, 
we can easily solve the original equation (l). 

As indicated previously, we may do an analogous procedure with 
system (6h), prodxicing a cascade of equations in the "opposite direc- 
tion”. Darhoux pointed out that if we followed the h suhstitution 
with a k suhstitution, we woi&ld not produce a new chain, hut would 
indeed revert to the original equation (l) . In fact, if we denote (l) 
hy E, denote the equations obtained from E hy use of the functions 
h, h^^, h 2 ^ 0.0 hy E^^, Eg# o.., and denote those obtained from E hy 
use of the functions k, k^, kg, ... hy E E_g, ••• ^ 

of equations appears as 

. . ., E_^, E_^, E, B^, Eg, . , . . 

If we shoiild take any E obtained hy use of an h , function, and 

n D.-1 

attenq)t to obtain a new equation using a k function, we would in fact 
produce equation E^_^. The same relationships hold for the E ^ equaticns. 

Darhoux discusses these and many others points regarding the 
Laplace cascade method, including the nature of the invariance of the h 
and k functions, periodicity of the h and k functions, and the form 
of the most general solution obtainable if the chain terminates after a 
finite number of iterations in either direction. It is the pvirpose of 
this thesis to extend this cascade method to larger classes of equaticns 
and to carry out similar investigations regarding these new applications 



of the method 






I 



6 



In particular we will first show how the cascade method can be 
applied to systems of n linear second order hyperbolic equations in 
two independent variables, with n dependent variables, discussing the 
invariant nature of our substitution functions, the general form of a 
solution when the chain terminates after a finite number of iterations, 
and the solution of some typical systems.. ¥e will consider the system 



( 10 ) 






n 

* r 

^=1 



a., 
10 



n 



3 ^ 

3 X 0=1 o 7 



n 



0 = 1 



®ij 0, i - 1,2, «.*,n , 



where a^^ ^ (x,y), b^^ = \j(x,y), c^^ = c^j(x,y) are continu- 

ously differentiable as often as necessary, and the 1 : 1 ^^ = Uj_ (x,y) 
are real functions of real variables. We will first express this equation 
in matrix form, letting A = ® ~ ^ ~ ^^io^ ' 



’"l" 


II 


'0' 


, 0 = 


00 • • • 0 


^2 




0 




00 • • 0 0 


• 




• 






• 




• 




• • • 


0 




« 








0 




0 

• 

• 

• 

0 

0 


L « J 


- 




L __ 



Then (10) may be written as 




+ 



AU^ . 




( 11 ) 



+ CU = (0) 




1 



7 



where U = 

xy 


^x ^y 
: 


<0 

II 


~>1 

^X 

« 




II 


57 

O 




o 

<y n 




o 

o 

n 






O 








a% _ 


i 




"TT 



Introducing the two D^rbonx n x n isatrix inTetriants 

H — A + BA = C , 

X J 

K — B 4- AB ° C 

y -> 

we will "be able to proeeed with developmemt of chains and solutionB 
by the cascade nsethodo 

Our seccsiLd extension of this nsethod will be to a single third- 
order linear hyperbolic equation in three indepesident variables, of 
the form 



( 12 ) 



u a(x,y, z )n ■+-b(xj,y, z )u +■ c(x,y, z )u 



sys 



ys 



Xz 






+• d(x,y, Z )u + ©(x,y,z )n f f(x,y,z )u +■ g(x,y, z ) u =. 

rt 0 

0 

Since we will be dealing in three independent variables, we will see 
that th3fee chains will resialt, in the ”x -direct ion”, the "y-direction", 
and the ”z -direct ion”. Also, sine© the equation is of third order, we 
will see that ei^teen "invariant” functions mast be introduced, ¥e 
will note that the invariance of these ftmetioas will nofc be of the same 
mtxsre as that of the corresponding ftmctions in the second order, two 



I 



I 






I 
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mxist be placed on the coefficients a(x 5 y>z), o * © # gCx^y^x) before 
a chain of equations can be developed© 

Finally we will generalize the results indicated in the preced- 
ing paragraph to the single linear hyperbolic equation of n^ order in 
n variable s^ which has the foim 






^X^()X2... ^ x^ 



n 






n-l 



u 



i = 1 



3 X]_ • © ^ ^i— 1 ^ 3Ci ^ Xjj 



n 



(13) 



ij 



isO = 1 

j 






n-2 



u 



a x^o . . dx^.i a x^^.1© . . ax . ^ . . ax^ 



o o o 



n _ 

" \ a u 

i = 1 



+ cu = 0 , 



In this section we will introduce an operator notation which will 

+Vi 

greatly simplify the calculations for n order equations© This 
will enable us to predict the number and form of the identities # and 
the number and form of the invariants corresponding to each identity 
related to the n^ order equation. From this we will be able to in- 
dicate the conditions necessary to produce a cascade of equations^ 
which will enable us to tell when a reduction in order will be possible. 

In our concluding section we will disciiss briefly the extensions 
made by Darboux, Dini^ and LeRoux to systems of second order equations 
in one dependent variable ^ to a single second order equation in n 
indepent variables, and to a single n^ 



order equation in two independent 



i 



9 



variables <. These extensions are treated more extensively in the works 
of these men noted in the bibliography. 



SECTION II 



AN ETTENSION OF THE PRODUCT INTEGRAL TO 
NON-HOMOGEENEOUS SYSTEMS OF TEE FIRST ORDER 



In Section III w© shall show how the cascade method, can he applied. 

to a system of second, order equations. The end result we hop© to obtain 

is a reduction to two first order systems^one of which will he homogeneois 

while the other is non-homogeneous , The solution of the homogeneous 

, 6 

system can he found by employing Volterra's product integral , As yet, 
however, this concept appears not to have been extended to non-homogeneous 
systems.'^ It will be necessary for us to do this now, to enable us to 
solve completely the second order systems. 

Consider the system of equations 



( 1 ) 



du^ 

dx 



^ "j + *1, 



^i(^) ^io ^ ” 1 ^ 2 ; 



where the a and the f . are given continuous, single -valued, bounded 

jX 

functions of the real variable x, on some non-empty Interval of the red. 
line, b ^ X ^ c, and the u^^^^ are n given constants. First let us 
write (l) in the notation of matrices. W© denote the row matrices 



U(x)~(Ui Ug . . . Uj^) 

F(x)— (fi f 2 . . . fjj) 

= (“lO “20 . . ."no> 



10 



11 



and the square matrix 



A = (^ij) 



n X n ' 



It will he necessary to assuzne that the matrix A is non-singular, 
that is that the equations of system (l) are linearly independent. 
In the notation above, system (l) can be written as 



( 2 ) 



dU(x) ~ 
djc 

U(b) — U 



U(x)A(x) -I- F(x) , 



Let be any partition of the interval /”b,^, such that 

^ point in the 

interval We then define by the relation 

(x^. X^.,)+ 

^ + ij + )(x^ - X^.^) . 



(3) 



Here I is the n x n identity matrix with ones on the diagonal and 
zeros elsewhere. From the form of this recurrence relation (3) we 
readily obtain that 






2rkil r m 



(M 






(x — X ) *l" 






Xm-1 ) « 











i 



3 



i 



12 



Next we eonsider a sequeiise of sueli partitions^ sueli that 

as m — ► ^ X ■= “ x^^i — * 0, Sinee all the ftmctione con- 

cerned are continuous, we laaj proceed to the limit; 



(5) u(o) =- “ "o I * + J F(^ ) [ A d-^ . 

■*“ ■ I 



Ih-»OQ 

e 



h 



h 



In this express icfli, I A= ^CA(^ )dJj1-l ) is the "right" product 

h h 

c 

integral of Yolterra, while is the row matrix of 

term. "by term Eiemann Integration ®of the elements of the row matrix, 

a, 

t 

In his discussion of product integration, Schlesinger proTed the 

8 

following identity; 



p- IM 



for p <. s < q,. 



Using this we see that 



c f c 

1--M 



h <• ^ < e. 



and hence 



( 6 ) 



r, 



Putting (6) into (5) w obtain 

(7) u(c) « u^ j A + ^ f(^)/Ja\ Ja 



1 



i 



( 



< 

I 



13 



Since 



r 



A is ccaistant with respect to the variable of integration, and 



since the scalar d ^ commutes with every square matrix, we may factor 
this product integral to the right, and write (T) as 



( 8 ) 



U(c) = 






_C 



This function U which we have derived is a fimction of the end point, 
c, of the interval If we vary this end point in any interval in 

which the functions a^j and f^ remain continuous, single-valued and ’ 
bounded, then U(c) hecames a function of a real variahle which we may 
call X, and hence 



(9) 



U(x) = 



^o + 





We assert that (9) is the solution of (2) and hence, element-wise, the 
solution of the nan-hcmogenous system (l) , To see this, we observe first 

— X Q /X\l 

that when x = h, T A = ■£ , while ^ F(^)/J A "5 = 0. 

Thus U(b) ■=■ U^. Next we note that each row of the matrix A is a 



solution of the matrix equation d U U(x)A(x) 



Therefore 



_d 

dx 




’(A). 



Thus differentiation of (9) yields 



I 



1 

( 



5 



I 



Ik 




— y(x) + u(x)a(x), 

which proves that (9) is the solutloa of (2). 

11 

Schlesinger showed that if A is a matrix of coostant 
elements^ then 




Therefore in the special case when (l) is a system with constant 
coefficients, (9) talces the more familiar form 




and thiB"BOlutlan can be veidfied even more readily by direct dif- 
ferentiation, We note the analogy between (12) and the solution of 
the single non«hamogeneous first order equation 



f 







I 



1 

I 

I 



I 




I 
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j- =■ a(x) u(x) ^ f(x) 



u(b)-= 



This equaticai has the familiar solution 



u(x) = 




_ C^a(’?)<3^ 

^ f(^) 47 ) , 



a form which we could obtain from (12) by taking the transpose. 



I 




I 



I 



I 

I 



I 

*1 



SECTION III 



SYSTEMS OF SECOND ORDER LINEAR 
HYPERBOLIC EQUATIONS 

A. Consider the following system of second order linear hyperbolic 
eq.uations ; 



( 1 ) 



X ^ n 

^ “i + y 









n 



J=1 ^ y 



n 

j =1 






1 > 2 ^ . . . n j 



where the coefficients hj^j, a n d Cj^j are all real functions of 

the real variables x and j, continuously differentiable in both 
variables as often as necessary, and Uj^(x,y) are real functions 

TP 

of the real variables x and y. To put this system in matrix notation , 



let A = (aij), B -(bij), C - (c^j), U’ 



we obtain 



V 



% 



(0) > r» 



and 



(2) 



xy 



ATI *1- 

X 



BU *h 

7 



CU 



( 0 ) 



i 6 











I 



I 





17 



where _ 



r 






i> 



m 

a xTy 



, etc. We will assume that the 



matrices A, B, and C are all non -singular. 

We introduce the two Darhoux matrix Invariants 







H »■ 




BA 


- C , 


(3) 




K — 


B u. ’ 

y ^ 


AB 


« C 0 


Then we 


note that (2) can he written in either of the forms 




'^x\ 




AU) 




4- AU'\ - HU *• 


(1^) 


^7 ^ 


. d ^ 


BU^ + 




^U 4“ BU) ~ KU — 


If we consider 


the suhstitutlons 










U 

1 




AU 


) 






^-1 = 


+ 


BU ^ then (4) becomes 






^1 

TT 


- 


HU 


- (0); 


(5) 




JL 

u_i + 


AU_j_ 


KU 


(0) ^ 



( 0 ) , 



( 0 ). 



If now either H — o or K— 0, where 0 



d 0 
o o 



e o 

o o 



our 



n X n 



I 

I 



I 



I 



I 



I 



I 

I 

I 

I 



18 



original system (2) will be reduced to two systems of first order 
equations, which can be solved by the product integrals of Section II. 
We will make the assiamption that none of our coefficient matrices are 
singular, except the matrix 0. 

Suppose that H S. 0 . Then we have the following systems to 

solve i 

+ BU = (0) , 

( 6 ) 

U = -t- AU . 

1 y 

The first equation of (6) has the solution 




where Y® (y) is a column matrix of arbitrary functions of y, and 

\ -1 

B( f ,y) d^ j is the inverse of the Vol terra product integral 
of B. Knowing U^, we may then proceed to integrate the second equation 
of (6), to obtain the solution 



(8) U(x,y) 



4 

where X° (x) 
If H ^ 




is a colimm matrix of arbitrary functions of x. 
0, but K = 0, we can solve the system 






AU 



_1 ' ( 0 ) ^ 



u 



-1 



u„ 



( 9 ) 



4- BU 






r- 



% 




«• 






i 



I 




I 

I 
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by the same methods utilized to solve system (6), to obtain the 
solution 



(10) 









Illustrative example No, 

Consider the system of three e(iuations in three unJoicwns 

^*^1 -I- y 4. y©^ ^ ^ y^i + y®^^i + j 

+ y ^U^. ^Ug - eyu^+ yug + eyu- =• 0 

^Jjy ^ 

^y^^Ul ^©^ ^Ug +yey ^Uj + ^113 +y^i ^-oyug + yeyu^ - 0 , 

^x^y ^x" ^"x Jy" 

This system can be written in the form (2) with matrix coefficients 

C A * 



A = 


1 


, B = 


H 

o 

o. 




y q7 




0 1 o 




y^ ^y ye^ 




o o 1 



If we compute the invariant H, we find H =0, Thus we may reduce 
to the two first order systems 






\ + 
X 

”7 + 



“l * (®) , 



AU = U, 



The first system of (b) has the immediate integral 

Ul = e"' 7(y) 

where Y (y) is a column matrix of three arbitrary functions of y» 
Using this expression for in the second system of (b), we obtain 



the solution 



I 

• roi^ 

• c 

f 

I 



♦ • • li 
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// \ -1 1 






— •*=* 



U(x,y) 

where X(x) is a colunsn matrix of arhitrary functloos of x. 
Bo Consider now successively the substitutions 



I 



U 

X 

y 



A u' j 

(x ) 



where A is an n X n matrix of twice continuously differentiable 
functions of x and y^ such that |AV^ o , and U 



-u- 



u: 



u^ = u^ (x,y). Under the change of variables U ® Au , ( 2 ) becomes, 

/y^x^y (A A y)u' ^ (B A t Ax)Uy t 

+ (CA A -t- B Ay ^xy) u'= (0). 
Multiplying from the left by A ^ 'we obtain 

^'xy + (A“\ A AV\y)u*^ + ( A“\a + A’^A )U« + 

X fl/ 



+ ( A'^c A + A“^a A^ A“^Ay +A“^ A (0)^ 

which is of the same form as (g)s 



(11) U’ + A'U^ + B'U'y + C’U* = (0)^ 







I 

I 

I 






\ 




m 










t 



1 
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(12) 



where A* = 

b’ » A’^BA^• A-^Ax^ 



A^cA+A"^aA +- A-^B Ay1“ a a. 






Let us ccmpute the Tulues of the carrespoudlng Darhoux iismrlaiits 
h' and k' for (11). 



H 



a’^ + b’a’- c' = 

( A ^ A A A*^^) + ( A A+ A A )(A aA'^'A Ay)- 

X X 



-( A "V A + A‘\ A + A"^ B A -J- A"^A •) ) . 

X y xy 



Consider 



AA‘^= I 



i,(AV^)= A^^ + ^A'^= ■)! _ ,0 

(13) 57 

- A‘^ MA'! 

• * ^ X 



Using this identity^3 we may differentiate the first term of h' 



to obtain 

s' ^ -/\^l\ A'^aA+ A'AAvA'^aA -A'‘A^A'Vy + 

+ A'^A^y 4- A'^baa A'-^ B Ay * A'^Ay A^A A + A Ajt A Ay ' 
'A'^Ayy - A^c A - A^A Ay - A'* B Ay - A'^A x„ = 

= A'YAy + BA-C)A 

= A" HA, 



In similar manner we may compute the value of 

k’ ^ b’ a'b'-c’= A'^kA. 

V 



I 

I 



i 










% 



I 

I I 







1 
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We see then that this change of Tarlahles U =* A merely produces 
a similarity transformation on the matrices H and K, and Indeed, If 
/\ Is chosen so that A S = H A > then H* *=• H, and should 
A K = kA > then K* =s. K. 

If we make the change of coordinates x y - 4^ (y'), 

then the resulting equation 

^x’y« •(- (ae)BUy« + 4^^(/)CU " (O) . 

also has the form of ( 2 .) » Here 



a' = iV'm a 

b' = B 

o' = (p'iJ) If' (y')C_, 

and the corresponding matrix Inmclants are 

h' ^ ^ (aO ip' (/) 

k’ (x) iP' (y) K, 



Finally the change of coordinates x « y’, y^ x’ merely has 

t ' 

the effect of interchanging the invariants, so that H * K, K =b H» 
Thus we see that the term '‘InTariant** is correctly chosen. 

C . In the likely event that neither H nor K is 0, let us consider 
the system 



(ih) 



Ui = 



u, -4- 

ix 




AU 






HU 

« 



Solving the second system of (Jh) for U, we obtain 
(15) , U - H“^ 1- h“V^ ^ 



Utilizing Identity (13) we further obtain 







‘P 
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(X6) 



-H 






U 



Ixy 



- H, 



+ E"^y Ui + h''^U. 



ly' 



Thea nnxtlipljring the first system of (l^^) from the left hy H; aad 
substituting (15) and (I6), collecting terms we have 



where 






A^Ui 



B U + 



= ( 0)^ 



•% 

(17) Bl - B , 

Cl = (HaS' - )B 4- By - H = A^Bi A-Bj^ . H . 



Thie equation ie again of the form of ( 2 ) and we may iterate our 
process in an attempt to reduce this system to two systems of first order* 
Computing the matrix invariants for system (17) we have 



®1 * *loc + Bi Ai -Cl 

•» (HAB"^ - HyB"^)^ + BCIAH"^ - ) - (HAB'^ - - 

«» By ^ H 

«2H-K«A3g4AB-BA + )B+ 

+ (HAH’\ - (Eyfl-1 )^j 



H. 

In a similar maimer we may solve the system 

w-l = 4- BU, 

B-OJ ■*■ AB.1 = KH^ 



2h 



to derive the system 

(19) = ( 0 ) 

where now 

A 1 - A 

“J- J 

B_3_ = EBK"^ 

C_i — (EBK"^ - K^“^)A - K + Aj. = 
“ ®-lA_i + ^ 

For (39) our matrix invariants will he 




K_J_ = SK - H - Bj, + BA -Jffi+A(KEK-1 - K^‘^) . (kbk'^. K k'^)a f 
+ (KEK-^)y -(jyc'^)y . 

If we designate the original system of equations ( 2 ) by E, 
we have been able to transform E into a system of similar form, which 
we shall designate E^^, by the change of variables % “ U AU, and 

j- gjr 

into another such system, to be designated by E_^, by the f'ha-n^a ©f 
variables U_]_ ^ U^; I^ovided that none of the resulting H or K 

invariants is non-singular, w© may ccntinue this procedure in both 
directions, establishing a chain of systems 

, . . E_g,E_^, E, ^ 

The question naturally arises, if after making the substitution 
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= Uy “h AU and discorerlng that H ^ 0 , we transform system, 
(2) into system (17)> what resialts if we then make the substitution 

? , Presumahly we shall again obtain a system of 
the form of ( 2 ), but what exactly are the coefficients of this new 
system in relation to the original syst«a ( 2 ) ? Utilizing the expres- 
sions for the coefficients of (19) and recalling that ve 

obtain 



^ 1-1 



“ HAH'^ - ^ 



~ HBE”^ 






^1-1 ®1-1'^1-1 ^Ij - ® 



= HCfj. + BA - BH“%y - AH”^ 4- 

1 IT -1_ tr“l ^”1 “1 ”1 

+ - h. 



Now consider the subBtltutlcn = Au*, and for A let us take E, 

l.e. ** BU', Tram. (12) we see the resulting system, again of the 

form (2), has the coefficl«its 



^.1 S + H-\ = A ^ 



b' * B, H + H"\ = B , 
1-1 X ^ 



i ' ^i-iS” 



Thus we have transfanaed into a system which is exactly equiyalent to sysSem 
( 2^, since 



i « 



♦ 
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E = l’ = E-^ 



HBE’^ = H 



1-1 



hence 



H *»H 



1-1 



This result may he Immediately generalized as follows; If, in system 
Ej we make the substitution we will obtain system 

make the substitution will obtain 

the system E^ for all i ■= ,,,, -2, -1, 0, 1, 2,,,. . The foregoing 
tacitly assumes, of course, that the corresponding Hp and Kp are neither 
e(iual to 0. 

Remark #1 ; We may define an exponential matrix where X is 

any n x n matrix, by the f omula 

^ 2 3 ^ i 

2 31 i»l 11 

How unless the matrix X cammutee with its x - derivative, we cannot 
X X 

assert that d * C dx , 
clx dx 



d^C I 4-1 (^Xf X^ 1 «^X4-X^(ax) 

^ 2^ax3ldxdx cLx * 



But if X=-Ay, say, where A is an n«n matrix of ecmstants, 

and y is a scalar matrix, then ^ * A ^ , hence 

ax ax 

X® * AyA^ ■= Ad^ « Ay ■ ~ X ^ L-( e ^) = , 

dx dx ax dx ^ dx dx 



Then consider the special case of system (6), when the matrices A and 
B are both matrices of constants. An integrating factor for the first 



system of (6) woixld be 



fedx ^Bx 
e'J = e . 



Then we may write 
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H- h'= 



-1 



== H, 



1-1 



hence H ** H 



1-1 



This result my he lumediately generalized as foUosrs: If, in, system 

Ej^ we make the suhetitution we will obtain system 

EjL^.;L, hut if we mke the suhstitution will obtain 

the system for all ±~ ,,,, -2, -1, 0, 1, 2,... , The foregoing 

ta^citly assumes, of course, that the correspcmding Ej^ and Kj_ are neither 
eq^ual to 0, 

Remark #1 ; We my define an e^onential mtrix where X is 

any n x n mtrix, by the formula 

^ 2 3 i 

© «I•»-X•^-X^-X +...-I+ X X . 

2 31 iTl ij 

How unless the mtrix X ccmmutss with its x - derivative, we cannot 
X X 

assert that d C * C dX , siiJce 
to dx 

d_C^=l4-i(dXXt X^ 4- 1 (aXj2 y ax X4- X^to) 
to 2 a? to 31 to to • 



But if X®Ay, say, where A is an n><n mtrix of constants, 

and y is a scalar mtrix, then to — A ^ , hence 

to to 

X^ * AyA^ = A^ « Ay r to X and thus d_( e 
tototo to^ to to* 

Then consider the special case of system (6), when the mtrices A and 

B ar© “both imtrices of constants* An integrating factor for tli© first 

system of (6) would be e . Then we may write 




L 
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Bx Bx 

e u, 4 e Eu. =■ (0) , 










- Y(y), 



vliere Y(y) is a colionn matrix of arbitrary fimstioas of y. 
Substituting e y (y) into the second equation of (6) gives 



~Bx 



Y(y) = Uy AU ^ 

An integrating factco' for this system is then ^ 

Integrating this system gives the solution 



0= 






Ay - Bx : 

C y(y)cLy -I- X(x)J 



-(Ay + Bx) 



y(y)ciy 

Ay 



[Je ^ ^ . 



where X(x) is a column matrix of arbitrary functions of x» 

Now consider system where no longer is H ^ O , and 

assume once more that A and B are constant matrices. We may inte- 
grate the first system of (l4) in the above manner to obtain 



U =. 





U^dy 



Z(x) 



>] 






Substituting this expression for U in the second system of ( 1^ ) 
leads to the equation 



U, + BU, * H 
■Sc ^ 







+ X(x) 



'1 



# 



Multiplying fron. the left by 




, this becomes 



Ay -1 
e H 




Uj^dy X(x) 
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If W0 now differentiate "botli sides witli respect to y we obtain 
( AH“^ - (U^ +■ BU3_) + 

\ jsy ^ 



Multiplying frcm the left by H e. 
finally obtain 



-Ay 



, ^and collecting terms we 



'’lx <-^1% <■ «1 "l = (°) 

where the coefficients A^, , and C]_are exactly as obtained 

for (17) . 

The curious thing about this second method for* starting a 

chain of eq.uation8 is that even if A and B do not ccmmute with 

their derivatives, we may still use the integrating factor and 

d <- 

operating formally as above, pretending that the equation ~ <S. .A 

were actually true, we will arrive once more at system (17) with the 
correct coefficients Aj^ , ^ snd « Thus we obtain a valid 

equation frcm an invalid operation, but an operation w'hich appears 
valid on the surface, and is completely analogous to the operations per- 
formed on a single equation of this type. 

D. If the chain has been continued to the i + 1®^ system of equations, 
we can readily establish the following recurrence relationahips for the 
matrix ^invariants i 

(20) %+l“^i -^1®! - ®iAi + Bi(3iAj^Hi - ) “ 



-1 



-1 



(^i'^i^i “ ^iy^i )B + (H )j£ “ (HiyE )x ^ 



H-1 i ) 



•*•2^ “X^ 0^ X^ 2^ • • • 
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These may he solTed for and to give 



(21) ^i+lj 



-1 - -1 






-1 



"l^ 



-1. 



i « e.c, -2, -1, 0, 1, 2,, 

Fran, these w immediately obtain the relaticias 



( 22 ) 



(23) 



-(E^A^% " )Ei +(HiAiH^ )^ “ (H;^^!^ )^ j 



-1 



-1, 



-(HjAjHj - HjyB, - (MjE/ _ 






0*0 “ “ Ju ^ 0 ^ 2^*00 

Finally (5) may nw b© -writton more generally as 



(21^) 



% + ’'l.l ^ 



Considering the first system, of ( 2h), suppose B is a matrix 
of constant elements o Multiplying both sides by the exponential matrix 
0^^ (see remark 1.) mid solving for w© obtain 



”Bx 



3x 



^i-1- ^i-l © ^ (<3 Ui) * 

o ^ 

This Immediately leads to an expression for U in terms of Uj^ and its 
derivatives with respect to x, namely 

-1 ^“Bx \ . Bx -1 -Bx N Bx -1 -By n n Bx 

E=H e L (e El e i_(e \ {...2£© %))..), 



( 25 ) 
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tedx vBdx; 

In the ease where the matrix B is such that ( © ^ ) = B©'* 

dx 

does not hold^the solution for U in terms of 'becomes much • 
more involved. Solving the first equation of (2l|-) for 'w© 

tor. %.1 = +B0j). Then %-2 = Ii.a'V-lj, + ®’l.l) = 

Hj. 2 ^ (El-l(El 3 |.+ BUj) •>- • (Oj^ + BU^)) . Iterating this 

process we see that eventually we will obtain in terms of 

Uj^ and the i invariants H, Using the relation 

^i “ ^i-l^i“lj 

we may, in a similar manner, obtain an expression for Uj^, in terms of U. 
E, We shall now give consideration to what may occur in the form of the 
invariants as we iterate along the chain of systems. Perhaps the most 
natural point to investigate is that of ’’periodic” systems, i.e. systems 
such that after J iterations, we obtain Ej =E. We shall say that such 
systems have period J. Consider first a system of period 1, Then 
®1“ % “ H, and »» K. Hence from (20) we have H * K, and 

H=.2H-E-A^ ^-AB-BA + B(HAH“^ - - HyH“^ )B+ 

f' 

In the special case where AB =■ BA and B(HAH“^ - HyH"^)-(HAH“l-HyH“^)B, 
this reduces to 

(HAH-l)i - - 0 ^ 

Which has the immediate integral 



g 

G| 

i 

I 

I 

f 



I 

I 

I] 



I 

I 
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(26) HAH-l -HyB’l - Aaa-Y(y)^ 

where Y(y) is an n x n matrix of arhitrary functions of y. 

If in addition HA « AH, (26) is further reduced to’ 

^ y(y) 

Hy * y(y) H , 

This may he further integrated hy the product integral to give 




where X(x) is an n i n matrix of arhitrary fxmctions of x. If we 
select for Y(y) a matrix of constant elements, D this solution 
can he written 

(28) H = X(x) , 

Suppose we have a system of period one and let us make a change 
of variables U » ^ u'. This will not change the periodicity or the 

equality of the invariants, for, as we have seen, if H * K, then 

h’« K /V = 

Let Us select A however, so that our coefficient a' 0 , This 

simply requires that A ^ = -A A , and hence 

A.), 

where X (x) is a column matrix of arhitrary functions of Xo 






1 

s 
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» I . » » ' • ' t ' » ' » ' • f > # t ■ t 

Then H=A 4-BA -C=-C«K, hutS + A B - C . C, 

I t ' 

hence ® . If we make further assumptions regarding the 

character of A and B, our system can he reduced further. Let A and 
B he matrices of constant elements. Then (29) can he expressed as 

•Ay 



(30) 



A = fs»^ * / \ 

f \ - O X (x). 



Since B*" ^B A + A"^At > also atten 5 >t. to find A such 



that B » 0 • This gives 



(31) 






e Y 



(y) 



Comparing (30) and (31) and noting that if and only if 

PQ « QP, we may take for our arbitrary matrices X* » 6*^- , Y* — ©"^y 
to get the matrix 



(32) 



A=e-‘V 






if and only if AB * BA, 

Using (32) to change coordinates under U = our resulting coef- 
ficients are v 



A-* ~ 0 j b’ — 0 j 



C ’ = -H* 



and our system (ll) has the reduced form 



(33) U* ^ hV_^ 

xy 

which is the matrix analog to the telegraph equations, Since^ if A and 
B are constants such that BA •»» AB, it follows that H =. K, we have 



proved the following; 










i 

’^1 



t: 






.1 



j 



\ 






) 
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Theorem I« A system of oquatlans of the form (2) havijig constant 

1 ^ - > ii-g^ 1 — w—^tm\m\mmm-\n-mmmmmm wiiiia m \m\ nmmammm — r — m—r 

matrix coefficients A and B can he reduced to the form 'D^ ~ H'U * 
hy a change of Tarlahles 13 ■s /\u' if and only if AB = EA. 

Illustrative example 

^ ^+.10U3_+45U2+Ilt9h3* 0^ 

dx ^x ^F" Tj ^ y 

2>2 ^^^3 = 2^Uj^ ^ll^Ug ^.:^2^U3 +20u3_ - Ijug -»-49u2 - 0 

-ST 7T 7T^ ^ TT 

^^If ^>2i 17^f 6oj^^ 19^12+T9U3=0 , 

W 5T 'hj 



This system may he -written in the form ( 2), with the constant 
matrix coefficients 



A ~ 


~1 


3 


5' 




Ba 


7 


2 


28~; 


C =’ 


"10 


45 


149 ” 




2 


-2 


1 








11 


12 




20 


-11 


49 




0 


1 


4 






2 


2 


17 


1 


6 


19 


79 . 




— 




- 


1 




L. 






] 


— 






We observe 


first that 


det 


A ■= 


-23, 


det B 


■=* 529 , and det C 


= - 8 , 792 , so 



that all matrices concerned are non-singular* (notice that B ^ A^ and 
C =:. A^ - I,) Computing the values of H and K we obtain 



I 



i 

1 



r 



3h 



H = +• BA 



C - 



0 0 0 
0 0 0 
0 0 0 



11 




lJf9 




'lO 


1^5 


1J^9 


20 


-16 


h 9 


— 


20 -17 


49 


6 


19 


80 


1 


6 


19 


79 




=■ I 



K=By'|-AB-C - 


'o 0 o' 




11 45 149I 


'10 45 149' 




0 0 0 




20 -16 49 1 “ 


20 -17 49 




0 0 0 




6 19 80 1 


6 19 79 




~1 0 o“ 


— 


I 

A 






0 10 










0 0 1 









Thus H = K = I, and w© further ohsorv© that AB » BA. The conditicais 
of theorem I Being satisfied, (c) reduces to 

(4) 



U. 






H U. 



Under the transformation U */\U*, however, w© know that 

h' ■= A H/\— /'S BO that (d) is in reality a 

system of three "uncoupled" tolo(.,i‘aph ©quationa, 

(.) = u’. 

Each equation of this system has the same solution u', which can 

15 



he found in the literature. 



Thus (©) has the solution U = Fu* 



t 

u, 

u 



Having detenained U , w© must then compute the solution U, using 
the relation U — A Equation (32) tolls us that 



A=e ' = 



-y 



e 



13 5 
2-2 1 
0 1 h 



- X 



7 2 28 

-2 11 12 

2 2 17 
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hence our solution is 

U 



-(Ax + By) 

e 



Un 



F« We consider next the systems of period two, so that E2 — E, 
and hence % « H and Kg — K. Equations (2p) and (22) tell 
us that 



(3^^) 



2H 2K - A13 . +A^Bi - B^Ai -J- Bi(KA^K"^ - K ^‘^) 



(KAiK’^ - KyK"^)B3^ 



-lx 



.-l^ 



and 



2K»2H-A3-4-AB-BA + B(EAH"^ - " 



(35) 



,-l „ „-l 



- (HAH"-^ - HyH"^)B + (HAH'^)^- (HyH"-^) 



r-l^ 



Summing these two equations we obtain 



( 3^ ) c — 1 ^ r 1 

+JBi KA^K - KyK ) - (KA^K - KyK )Bit + j B(HAH"-^ - HyH -^) 

-(EAS'^; - B,nf)B^ ' - {E^f)j = O , 

Looking at the terms within each pair of curly brackets, we note that 
under certain obvious conditions of commutativity this system reduces to 

(37) . 

Equation (37) IIS'S "the iranodiate integral 

KyK"^ Eyf^ = Y(y) 



(38) 
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where Y(y) is a s^iaare laatrix of functions of y. 

If wo ln5)Oso another condition of caEnmtativity, whereby 
EH ■* HE, E^ SSL. EKy then an integrating factor will he multiplica- 
tion frcm the ri^t hy EH. For this loads to 

EyH 4- HyK - Y(HK) ^ 

J. (HE) =■ Y(HK) . ; 

Hence 




where l(x) is a square matrix of functions of x. 

These functions of x and y are in fact determined hy the coeffi- 
cients A, B, and C. Thus we have proved Thecn^em II ; 

If a system of equations of the form (2) has period two, and 
the following conditions of eummutativity are satisfied ; 

(a) E/^ =. A^E 

(h) HA = AH 

(c) AB = BA 

^ ^ 1®1 “ ® 1^1 

(e) BiEyE“^ •= K^’^1 

(f) BHyH"l = HjH"^ 

(g) EH = BE 

(h) E^ - HE^ 

then the product of the two matrix invariants H and E has the form 
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HK » Y*(y) z’(x)^ 

where Y* (y) and x' (x) are squaro matrices of fimetions of y and 
X respectively. 

Since H is non-singular, there exists a matrix e defined hy 
the relation H « 0® , so that < 3 > = log Now if the conditions 

of theorem II are satisfied; and in addition, (39) takes the simplified 
form EK ~ I, then (34) can he written 

SB =,2 K - , 



or 



hence 



2 H - 2 B"^ = (B"^ ) 

y 



m 



2 £. 



-ae. 






Furthermore, if O has the property that © 
then B ^ — E , so that 

^y Ty 

(^1) ^ g-lg 

77 y 

Substituting (4l) into (40) yields the eq,uation 



^0 ^ 
Ty 



(1^2) 



^ ^ - 
^3:^y 



(0 _<g> 

= 2e -2e. f 



and this non-linear system of eq.uations will then have the solution 

© = log E. Eq.uation (42) is the analog of Llcwville’s 
17 

egiiation. ' 




I 
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G. Wo ccmo now to a conaidoration of tho general form of the solution 
■when scaao % =0 0. . Tliat la, wo iterate the process 

forming the chain of equations E, Ep E2, . . . E^, until ss. 0 , It 

will then he inq>ossihlo to form equation E^ ^ hence E is the first 
system for which the H invar iant vanishes* The first equation of (4) 
then has the form 

A-b^l +■ *1 Bi / 3 Di + AjuA = (0) , 

Hs? } \i- > 



This has the immediate first integral 




Y(y)^ 



whore Y (y) is a column matrix of functions of y only. Upon inte- 
grating^ a second time, we obtain the solution 

U.U,y>*^ Aj^(x, ^ )d)"^ ^(x) f: , } 

whore l(x) is a column matrix of functions of x only. 

This solution is of the form 



jl(x) 



-I- 




J 



where and 0 are n x n matrices of certain well defined functions of 
X and y. 



Using the iteration scheme of paragraph D, wo see that the solu- 
tion for U has the: form 

(44) U* + Je rdr) + ( I + feidy) + 

^ Yd 7) 
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designat® certain well defined n x n matrices of 
fxmctions of x and y, and are eciaposod of certain products of 

B, and their derivatives. We note that the coltaan matrix 
Y always appears mdor the integral siga, in the most general case. 

If the houndary conditions are such that Y ®(<i) we see that the 
solution is of the fom 

® 

and this is the moat general solution in which appear arhitrary functions 
of X, and no integral signs. 

Conversely, if the original system of eq.uations has a particular 
solution of the form of then the repeated application of the cas- 

cade method will certainly load, after a number of iterations not greater 
than i, to a system for which Ej ^ ^ , We will now prove this, i,e, 
if (45) is a particular solution, the number of iterations necessary 
before any H-invariant vanishes j is not greater than 1, If we substitute 
(45) into the system (.2), we will obtain an ei^ression of the form 

yft-i Vii' <-...■(• (o'). 

From the arbitrary character of X, we assert that 



(46) ®0 j — 0,1,2,,,®, i+l ^ 

If we compute the form of and obtain 



^ ^ y y 

>14 1= A<2 

a y 






o x^ y ^ ^ Xy ^ ^ 



Eq.uation (46) ingilies that 



(^T) ^CLl ^ . Q. - 0 

-jj A 0 ^ 

thus - ^^i-1 ^^1-1 >1- ^ 4- ^ 

^ 7 ^ 31 \3y 

-A+BA-cAa^ - -1^1=0, 

^ “T7 

Eance ve hare that 

(48) ^i.l -I- . 

ly 

Recall now the first substltutioa 



(49) 




AU 

4 



If we substitute the expressicsn for U given hj (45) into (49), 

sJcd Jrtlliz©, (?t7>), we see irmediately that the resulting eaqpression for 

Ul will have no derivatives of order greater than i - 1. In addition, 

(i-l) 

(48) informs us that the coefficient of X ' in the expression for Uj^ 
will he Whenever H vanishes, therefore, the coefficient of 

in this expression vanishes also. That is to say, i.f H vanishes 
the expression fcxr would have no derivatives of X of order greater 
than i-2. As a conseq.ueneo, repeated application may load to a system 
Ej, for J < i, for which the invariant Hj = O j otherwise we will 
obtain an equation Ej^, which has a particular solution of the form 

Ui = 0. X . 

®o swnmarizo, wo have proved the theorem III ; 




I 




^ ^ than the mathod of Laplaco vill leads after a ntaaber of 

ItoratloiiB not greater than i;, to an oguatlcei yhich is lirfcegrable . 

H. Considor nov ©xpa:“9ssicai0 of the fom 

I (1) 

U a 0. X •«- + ... f a^x' J 

which contain a matrix of arhitrary functions of x^ and deriTativos 
of this matrix up to a specified <hr>der. It is obviously quite often 
possible to express U in terms of derivatives of an order neater than 
specified. That is^, if X is expressible as a sum of matrices! 

where (8j ^ ar» certain matrix functians of x, and X^ is a 

matrix of arbitrary functions of x/ then U may contain derivatives 
of Ii up to order i . Conversely it may be possible to reduce the 
order of the highest derivative appearing in U. For instance, if 

U « (X 4 X^) -J- 6(x> 

then the substitution X^ ~ X 4- x' will reduce the order of the 
highest derivative by 1, 

We will say that the matrix U has rank i41 with respect to 
I if the highest derivative of X appearing in the expression for U is 
of order i, and it is iorpossible to reduce this order by substitutions 
of the type described above. 

We assert, that if the system (Ei) is the first for which the 



inTariant vanishes, then there exists a peirticular solution of 

the original system, con 5 )osed of a column, matrix Z of arbitrary func- 
tions of X, and derivatives of this matrix up to end including order i. 
This solution is irreducible in the order of derivatives, and hence has 
the rank i 4’ !• 

For if there were a substitution which would reduce the order 
of the highest derivatives to 1 -/^ , say, then there would be a 
system (Ej ) J = i < 1, for which this new expression would be a 
particular integral. But this woul(^ In^ly that Hj — 0, which is con- 
treury to our hypothesis. As a result, it is evident that if (E^) is 
the first system for which 0, then no other system (E^ _ j^) with 

positive or negative indices will admit a solution of rank i+1 with 
respect to i. 

The results of the preceding discussion apply without modifica- 
tion to the second substitution, l,e, the substitution which results in . 
the chain (E_]^) (E_ 2 ) , , , If this chain eventually terminates at some 
system (E_.) for which K . — 0, then there exists a particular integral 

d ""J 

of rank j + 1 with respect to y^ 

(J) 



u*- 6y + + ,,.,+ej Y 



and conversely. 



I, We are now in a position to construct aH the systems of a given 
dlmenslon^where by the dimension is meant the number of dependent vari- 
ables)of the form (2), which will lead to a general integral by this 
extension of Laplace's cascade method. Suppose, for example, we wish to 



construct a system of dimension n, which will have a solution of 
rank i 4 -l with respect to x« Then the chain must terminate after 
i operations, so that the system (Ej^) will he integrable. First 



we choose eurhitrarily matrices and non -singular, of order n. 

Then the eq,uation. 

(50) % =• ^ B^Ai - Ci ( O ) 

will determine C Then a particular integral is given hy ( 43 ) o 
The value of the invariant for (Ej^) is determined hy 

( 51 ) % ~ -t 

and then the relations (17) and (21) will permit us to calculate 
the invariants for the systems (Ej^ ^3^), . . 0, (E), 

Viz; Hi„i * 

®i-l “®ij T 

-1 -h 

^i-l = (Ai (4- Bi-ly' 

(52) = Hi.3_ Aj Hi.i + 

y y 

These relations also permit the calculation of the coefficients 
A, B, and C. From the expression in ( 52 ) for and Bj^_3^# we 



can show that 



Mj- 



A » (H H . . H^H) + 

^ (it 



B = B. 



and, having determiiied K, 

C a By + A3 - K , 

Having thus determined all of the Invariant matrices H., Oi:J ^ 1, 

J 

We may then compute the general integral for U hy the Iteration 
process of paragraph D. 

We may determine in a similar manner all systems which ter- 
minate in “both senses, and hence admit a particular solution of the 
form 



which contains no integral sign. The pMceding expression has rank 
i-l-1 with respect to x and rank J -tl with respect to y. The sum 
i + J is called the characteristic nujiber of the equation. This 
number does not change upon successive applications of the Laplace 
Method. In passing from system (B) to system ) for example, 

the ntimber i is diminished hy h, hut the number J is Increased 
hy the same amount j the sum of the two is michanged. This is apparent 
if we consider the system (E^), with invariant Ej^*0. This is of 



the form 




i 

IL 



/ 






which admits a solution of rai)k 1 with respect to x, and of rank 
i-hj -H with respect to y. For "brevity, let n » i J and con- 
sider the substitution 



where 0 is a colunai matrix of unknown functions. 

i® - Ai ♦ ( lAidy)"^® 
d y 3y 

V 7 ^ 

Using this (53) becomes 

(54) j + B, J (jA^ay)-" ^ = <°> • 



An integrating factor for 
JSj^dx. Thus we may write 



(54) is multiplication from the left by 



( 55 ) 

Denoting 



jBj^dx • (|Ady) 

A [oc 
I 



(lAidy)-^ = 

— 1 *1 

“ ^ ) (55) becomes 

^©1 = (0) 

TFi 



( 0 ) 









whicli may "be Immediately integrated to give 

( 56 ) ^ = e = ^ 

where > designates a column matrix of arbitrary functions of 

y, and X 2 a colximn matrix of arbitrary functions of x. We know, 
alrea( 3 y, that ( 56 ) admits a solution of the form^ 

( 57 ) 0 - I 1 + 6 Y i- 6iT^ + 

in which are certain v^ell -defined matrices of functions 

of X and j, and Y is a colm matrix of functions of y oixly. 

If we substitute ( 57 ) ( 5 ^) and take an arbitrary numerical 

value X, we find that Y and Yl are related by an expression of 

the form 



(58) Y^ n + 1 ^ 

in which A , A ... A certain sq,uare matrices of func- 

tions of y only. Substituting the expressions for ® and Y^ given 
by (57) and ( 58 ) into the first equation of ( 56 ) implies that 

(59) i.[6l+(BiYV...<- ^yH = OC ( AY + A3^rV...f 

In order that (59) be a valid equation, the coefficients of like order 



1^7 



derivatives on each side must he e(iual. ,In this manner we obtain 
the system 



( 60 ) 



^ = ocA 

Jy 

|6.i + e = <^A 

d7 

u 



1. 

+ ®h-l “ ^ j 



6n - ^n*l . 



(61) 



If we eliminate the from (60), we obtain the equation 



a linear equation of order nfl with respect to .OC -who^ matrix 
coefficients contain functions of y only. If we solve (6o) for 
the we obtain the values 

(62) ®’n=0CA„*i^ 




which permits the determination of the expression for O • 

The relation (58) , between the matrices Y and Y^ ) per- 
mits an arbitrary clioic© of either one or the other for the defini- 
tion of the value of © , by either (57) or the second equation 
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of (56). The choice of the former, however, will give the most general 
form of the integral for @ , 

As an application of the preceding discussion, which is com- 
plementary to our original problem, we see that it is sufficient to 
choose n + 2 square matrices of arbitrary functione of y, namely 

^ solve the resulting linear equation of 
n -I- 1^ order in one ind-ependent variable, (61), for the matrix OC , 
in order to determine the form of the equation (Ej^) * The formulae 
(62) and (57) then provide the general integral of the equation (56) 
and hence ( 53 ) • Repeated application of the substitutions of Laplace 
will enable us to determine the system (E) which has this general 
Integral. 

How consider an arbitrary linear equation of (n 1 ) order of 
the form 



/ ^ (nfl) 

( 63 ) MO( + M-OC+...+M OC 0 , 

1 n + 1 ^ 

where M, are any n+2 square matrices of known functions 

y and 0^ .is a, square matrJjc ©T functions of y . l.et A, ... ^ri+1 
be defined by the identity 

(n^l) 

A.W +• /\,W + /\,^ W 4. ... 4. /\ V 

( 64 ) 



n^l 



W 



.i. 1* 



Aw +• Aj^w + A2 w* * + ... A 

, .n+l ]\n 4 -l 

5y ‘ ^ '"2"' "•••+•(- 1 ) ^ 

matr 

then determines the A'® terms of the M's, f^ example 



-MW - _^(MjW)+ ^ (MpW) -... + (-l)'''^^ 

Tr 3 ^. , 

where ¥ is any arbitrary square matrix of functions of y. ^ 04 ) 







1 




e 




I 

I 



k9 



,=M. 



( 65 ) 



^ 7 






\*1' t-i) ^M 



3 y 



n-t-1 
n+1 ^ 



ml 



Iix llgib.'b of ©(luatious (60) and (65) ‘w© may se1 



M 



n-j- 1 



(- 1 ) 



n4- 1 



0? 






That is, if OC is a solutiou. of (63), then w© may determin© the 

fxinctions Aj by the identity-* ( 64 ), and the "by th© eyst©m. 

d 

(60). Upon substitution into (57)# w© s©e tliat th© general integral 
contains the column laatrix Y and its derivatiyes up to order n. 
Thus w© have proved th© following; 

Theorem IV; In order that the linear eq.uation 




admit a solution of rank: n4-l with respect to j, that is for th© 
general solution to he of the form 

O- X 4- Y 40, y'4 ,.o 

i. O u 

such that the solution cannot he wltten in analogous form in which there 
appear fever derivatives <gf the arbitrary column matrix Y , it is 
necessary and sixfficlent that OC , considered as a function of y^satisfy 
a linear equation of oic'der n 4- 1 whose coeffidents are functions of y, 
and that <3C does not satisfy a similar linear eguatlon. of lesser order* 



Jo Th© results given in the section have been kept in close analogy to 
those contained in Chapter II, volume 2 of Lecpns Sui* La Theorie Qenerade 
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Dee Surfaces by Gaaton Darboux. In this chapter Darboux acconqjliehes 
more with the single secoed order linear hyperbolic eq,iiation, than is 
available to ns with a system of such eqnaticaas, due to such causes 
as the nan-commutativity of the ring of n x n laatrices over the field 
of functions of two real variables. It is of course probable that some 
of the results obtained by Darboux which were not considered by us in 
this sectios.^ could be obtained far systems by matrix methods. The 
reader is referred to this vezy fine work by Darboux which has been the 
inspiration for much of this thesis. 



lY 



HTtT5FJ]0LIC EQUATIOIfS 0? 'IHIRD OHDKE HI 
riiOT BrosLtSKBS!:'’? 'muMiMs 

A, Let us now tiira our attontisn. to the iinoa'*” oqLmtl.onB of third 
order in three .indeptuidenfe Tay,'ia'£iles vhlch hare the form 



(1) (?f(v) 



u, 



X 7 ii 



siif U* 



+ ou 



xz 



C’ 






+ 






f a 

z 



gu s 0 



The coeffici.entB L, d,„ e_, and s' are tc he considered as 
functioaa of 2 :^ and c oxit Lt-uoae d'.!f'feroiTi:tiaS>j.o as many times as 



wo maj nsedo Wa wieh zo atoacs' thJ=3 ;^'r>:-'h;..eat .lr.i the max'ner of Laplace dis- 
cussed preTiousi^, Axi the h-xpe .rf '•■oduclAg (l) to a sjsteai of three 
first order equations „ failAiig ir ~nis ‘‘Xs will then attempt to cascade 
the equations^ jhi tho hope tlat afto:" a f number of iterations the 

chain will t©2"mJj:jato with yt>nd.s.V.ir£ Iwerlajatb^ and that the i*©sulting 
system can then he reduced a» vas original?^ desiredo 



To coDiSjenee this opox’aticc'i^ w« imist considei' puhstitutions of throe 
different typos, naToly 



(2) 


XX. - 

X 






au^ 




(3) 


xC - 

X 


V 


•r 


TbU-, 

1; 




(4) 




ys 




h-a 4- 
z 


cu 4- du« 


Let 


n& Btax^t with ot 


' rar 


’isoles ("), 


Uaing the definition of 


the 


operator oC which ap.w’ixre 


in 


(1). 


'ifiri I'oad. 


verify that 
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(5) 



u 4- 4* eu]^ + ^ (u) 4" + 

^ Z 2a y 




where 

( 6 ) 




a 4* ac - e 
z i 



a 4* ah “ f 
y 



4- ca + ha +■ a 
y z yz 



S 



are defijoed as the first three x-lamrian.ta. It is apparent fi'cm the 
symmetry of the equation that we eould have Just as well made the substi- 
tutions 





( 2 ’) 


^2 


— 


“y ^ 


bu 


or 






( 2 ”) 


"3 


rr 




cu. 






These 


lead to equations 


similar to 


(5) with the 


y-lnvariants (arising 


from 


(2')) being 


hg 


— 




be 












— 


h 4- 

X 


ba - 












— 


h + 
xz 


cb + 

X 


ab 

z 


4- be - g^ 


and the z -invariants (arising from 


( 2 ' ')) being 












3 ; 


"y 


cb 










^3 


— 


°x ^ 


ca ~ 










m 

3 




C 4- 
xy 


be 4- 
X 


acy 


4* cf - g 

• 



It is also apparent that whatever results we obtain from substitution (2) 







I 
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would "be symmetrically obtained ft*cm aubstltutioaiB (£') and (2**)# 
hence at this point vte will limit the majority of the diecusslon to the 
effects of substitution (2), 

Should we he 3?ewarded hy finding = h^^ = = 0, then (5) 

will become a second order hyperbolic equation In only two TOrlablos 



We may then apply the Laplace method to (7)/ us described briefly in 

18 

Section I; or in more detail in the works of Darbooz. 

Consider next the change of variables (3)^ which leads to the 
equation 



symmetry of the equation indicates that the foUovlng changes of variables 
would lead to equations similar to (8); 



(7) 



u, bu- + Ctt- + du- » 0 

■yz z y 









(3’) 





i “ "3, ’““3 * 



These substitutions lead to the respective l3)vaxlants 








sz 


a 

yz 


-{- (ac) 

y 


H- 


ha 

z 


+ ahe 


' 


■4 


ss. 


xz 


+• (ah) 

z 


+ 


c“b 

X 


+■ ahe 




m3 

2 


— 


1) 

X2 


(he) 

X 


4 


ab 

z 


4 ahe 


“ s 


ml 

m 

3 




C 

xy 


■<” (ac) 

y 


4" 


*bc 

X 


4 ahe 


- g 

) 


2 

m 

3 


— 


c 

xy 


4 (he) 

X 


4 


ae 

y 


4 ahe 


- g 

0 


In the case of equation 


(8), 




— 


^2 = 


II 


0, we will 



have reduced our system iimaediatoly to three first order eq.uations 





^X 


+ 


au — 


"lo 


(9) 


\ 


4 


hUj^ “ 


2 




•^Z 


+ 


* 





which can he solved in Inverse order hy quadratures. 

Thirdly we consider the change of variahles (!|-), which leads to 
the equation 



( 10 ) 




au_2. 



0^ i^) + 




-I- k^u^ -I- u. 



where m_]_ = 4- ad - g is the first - x -invariant . Once again 

ayinmetry consideratians show that the change of variahles 



and 



^-2 = 




4 


au^ 


4 




4 


^-3 “ 


u 

2y 


4 


au 

y 


4 


hu 

X 


4 



eu 



j 
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lead to tile -y -invariant and -z-lnvarlant 



m_2 






eb 


- s. 


m.3 


- ^z 


4- 


fc 


- s. 



respectively,, 

If we are fortunate once more, stich that ^2 * ^ ^ ^-1 ~ 
we may solve (lO) hy quadratures, and then apply the Laplace method 
to the non-homogeneous equation (li-) . 

B. We have in the discussion of paragraph A encountered some eighteen 

% 

expressions which we have labeled "invariants". Organizing these in 
somewhat more orderly fashion, we have the following tables 



hg 




I’z 


+ 


be 


- -d 








^3 


— 




+ 


be 


- d 








k 


S 


a 


-1- 


ac 


- e 








1 




s. 














k. 


d 


c 


+ 


ac 


- e 








3 




X 
















r 


V 


+- 


ab 


- f 












b 


-1- 


ab 


- f 








*2 




X 














“l 




a 

yz 


+ 




+ Cy 




ad 


- g 






xz 


1- 


ab 

z 


+ cb 

X 


+■ 


be 


- g 


m3 


- 


xy 


+ 




'”=1 


f 


cf 


- g 



f 

1 

] 

D 
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= t ad 


- 


8 




“^-2 


so e +■ he 
7 


« 


8 




“^-3 


— f^ 4- cf 


“ 


8 




4 


- V + ^ 




4- ahe 


“ 8 


4 


“a + (ac) 4 
yz y 


ha 

z 


4" ahe 


“ 8 


4 


“ h 4- (ah) 4- 
xz z 


ch 

X 


+• ahe 


- 8 


4 


= h 4- (he) 4 

XZ X 

e 


ah 

z 


4- ahe 


- 8 


3 


« c 4- (ac) 4 
xy y 


he 

X 


4* ahe 


- 8 


™2 

m 

3 


^ c + (he) 4- 
xy X 


80 

y 


4“ ahe 


- 8 



W® now ■wish to investigate the character of these eighteen ** invariants", 
as was done in paragraph B of Section III for systems, in order to see 
if the t®na "invariant" is appropriate. That is, we are concerned with 
the change in these espressicns when we make the various changes of 
variables 



( 11 ) 

( 12 ) 



u 

X 

y 

z 



X 



u, where X ~ X(x,y, z) is at least 
twice continuously different iahlej 



= ^ (x’)j 

= 4 ^(/), 
= Z.(a')3 



X 



f 

7 ^ 



(13) 






f 










4 



I 
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Let us coneider first the change of variables (11). Eq.uation 
(l) is then transformed into a new equation of the same type 



ilk) <;^''(u’) - u^2 ^ ^'^'xy"*" ■*' ®'V^ 



+ f U2 f 



g u =0 



but with new coefficients 



(15) 



„ . 1 H 

^ + X 77 

“ b + ^ 

^ X 77 

= d+s.^fS.lh.<-i 

X d y Adz ^ dyd Z 

_ ^ , c j. a , 1 

X XTT > <^xVz 

2 

= f + ^ + s. 1, 1 ^ ^ 

X 77 X ^y X Txjy 



g'= a I l^d ib <^^A+C^^A . ^3X ) 



/ r \ 
V • ; 



(16) 



Using these new coefficients, let us ccmpute the three new x-invariant, 
k^, 1^, and . We find 

»e a^ -r a c - e ~ ^ 

ij = + a’t’ - f’ = 1^^ 

“i = a^2 ■•■ *>'.'2 * a'a^ + a''^' - s' = m-i + jX ^ ll j X 

X ^y X V ^ 

The new - x-invariant, becomes 

(IT) Kill “ - s' * “-1 + ^ ^ + 1:^ , 

X 7y X ^ 

o 2* 

while the new invariant, , 



becomes 
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( 18 ) 



2 “ 

“1 



^7Z 



+ (a'-b’)„ 4" c'a’ t 



a 



S 




Tss + yir ^ j-T^ 



Froea i^yiamstrj, similar resialta ara ©btaiMd for tli© other 
h*, k', 1*, and m' inmriantso Tims to to® that th® h^ and X 
invariants, whieh appear as eoeffieients of th® x, j, and z derivatives 
of u, r©spe©tlV0lj, are indeed trrj® invar ianta rader th© ehange of 
variahles (ll|. Th© m invariaasts, era, th® other hand, do not repro- 
due® themselves ©xaetlj mder this ehang© of variahles, hut instead re- 
produce themselves plus a linear eo^'ination of the true Invariants. 
Hence wo shall in th® futOT® refer to thee® m iavarianta as q.naei- 
Invarlants . (Those will appear again in Section T.) 

If we should nak® the change of variables (12) and compute the 
result!!^ invariants, we will find that 



k;[ a , 

\ “ j 

( 19 ) ^ 

Finally w® not® that (13) merelj inter&angee the x-invariants 
with the y -invariants^-, the y-invarlaats with the z-invariants, and the 
z -invariants with the x-invarianta . 



C«: Let us now develop the cascade of equations in the zcamssr of laplace 

and Darbooxo Consider first equation (5) th® event that at least qm of 
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the invariants is not zero. We wish, to transfcxrm (5) into an equation 
for Uj which is in the form (l) , (See in this connection remark in 
footnote (5)) In order to do so, we proceed to integrate equation (2), 
obtaining • 



( 20 ) 



u e 



-X 



adx 



f 1 

\e Uj^dx + X( 7 , z)l 



where X(y, z) is an arbitrary function of y and z only. Then differ- 



entiating 


( 20 ) with respect to y 
- \ adx ^ 1 


and 

■f 5 


z yields 
.adx 


(21) 


u = e 
y 




r 


U]^ dx + X 




- ^(iadx) 
^y 


e 


Jadx 


RK u, ax 



. \ 

(22) u_s a ^ & 






- \adx- ^adz ^ 



(adx 

.s'* Uj^ dx X 






(23) 



Substituting ( 20 ), (21), and ( 22 ) into (5) we obtain 

-Cadx f r tdx , , 

I j ® Uj^dx + X »■ 



"lyz 

-jadx 



k 0 






> 1 , . 



dx 4 x; f 



■I 



^adx 






Uite + ij - 4 l|j. 



Now ( 23 ) is an equation for u^ alone, but it is an integral- 
differential equation, interlaced with arbitrary functions and their partial 



6o 



derivatives, and hardly in a manageable state. Our aim is to reduce this 
to an equation for u^^ in the form of (l) and to this end we will dif- 
ferentiate (23) with respect to x. Before w® do this, however, it will 
be necessary to make certain assumptioEis regarding the invariants and the 
coefficients. If these assumptions are made, the integral terms and the 
arbitrary functions will disappear when we differentiate. We will con- 
sider a number of different methods. 



METHOD 1 



The first assumptions In this method are 

(a) 1^1= li = “1 # Oj 

(b) a is a function of (x) omly. 

Under these hypotheses, (23) takes the fcam ^^ij. du-^^ = 

jadx + X + je^^^u^dx ^ ij ^ ^Ke^^^Uidr 4 t- 

j r V. J 

\ adx 

Multiply both sides by , and differentiate with respect to i to 

“1 

get 



= mj^ e 






adx 



u, 4 bu, 4 Ctt, . 4* dui ) - 

■Vz 'V 



(2U) 



- (Hi 4 bui 4 eui 4 du, ) 4 

^x m^ "Tz 



Jadx 

I 

“1 



— a 



(“2^, + ’’“la. + 


cui^ 4 du^^ 4 


®x^4 


adx 


(adx I 


[adx 


u^ 4 


^ u- 4 e '■ 


^ u, 



6l 



- \adx 

Kow wo may moltiply ‘both sidos of { 2 k) "by m-j_o collect 

terms to obtain 



u, -(- (a 
xyz 



^log mi > 11 . 






4 (c(a - ^ log m^ ) -h - m^) Ui^ 4 t(a - ^log mi) 4 - mi)ui 4 

TT” ^ -“Tir 2 

4 4 (a - ^log mi) 4. d_ - m,) u-, — 0 , 

d X 

Designate the coefficients as follows: 

ai » a - ^ log mi ^ ^1 ® ^ i ‘^l * ° ^ '^l ~ ^ 

®1 - ® 1^1 + ®lx - ”^i ) f 1 = l^iai 4 ti^ - 

®1 ~ ^l®"! + ^Ij^ " %• 

Then ( 2 ^) is in the same form as (l): 

(25) 4 a^ui^^ 4 hiu^2 *•“ '=1%^+ 4iU3^+®i;«ly+fiU3:.4SiUi 4 - 0 

Our now x-invariants are, frcm the substitution u^ =■ Ui^ 4 ^1^1> 

^ - aiCi ai^ “ ®i = ^1®1 ^1^ “ ^1°1 ” ®li + % 

= ®lz ” t mi s 
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Also 



1 , s + 



a 



V 



S a^bi- ^^log 

*** ^X*^X ^ ^X ” *^x~ ^ 



Ir (^x + ^1 ) , 

^ f a. 



y lyz 



gi 



‘JzJT’ 

_ ^ ^ log lUj^ 

^x^y ^ z 



=1 



og 



— m - ( d + b ^ ^log ^ c ^ ^log ^ ^ log m. 

“ — ' — ' — — - ^ac^KSE^^C=JttMK» T ci_ji^_»_ jj._ , .1 I , 

^ Z^ X 



If now the new X“in'mriants are all zero, we may proceed to reduce 
the e(iuation as described in paragraph A and solve. If, however, at 
least one of the invariants is not zero, we are in a dilemma, for the 
hypothesis with which we transformed (5) into (25) are no longer valid in 
the general case. We can no longer continue our chain. To enable us to do 
so, we must change our original hypothesis to the following; 

(a) * m^ ^ 0^ 

(b) All coefficients are functiongof x only^ 

(c) b ^ c S. d ^ 

Under these hypotheses, our invariants are also functicms of x only. We 



I 










y- 



I I 




r' 






« 
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may transform (5) into (25) as tefore, and the coefficients will have 
the same form. The x-invariants will he different, however, and in fact 
will he 



ki^ = ki . c^ = ki - d^ ^ 

(26) U * ll ■ \ " ll ■ 

Hence kjj^ = = m]^ , our hypothesis are again satisfied, and the chain 

may he continued. 

Observe that, as our chain progresses, o\ir x-invariants will always 
have the form 

V ' J =0,1,2... , 

Thus, if our invariants are not originally all; lerp^ they can only hecome 
zero if they are an integer multiple of the partial derivative of d with 
respect to x. 

Hote also, that to have a chain of y-invariants, or z-invariants, 
our hypotheses (a), (h), and (c) must he modified accordingly. 

METHOD 2 

In this method we assume first that 
(a) = 0 , 

Under this hypothesis, equation (23) takes the form 








% 




6k 



(27) '*■ 






i-adx 



Multiply "both sldee "by e 









adx 



Uj^dx +• 



■]. 



adx 



and differentiate with respect to 



mi 



X to get 
^adx 

Q (fll “ ^ log OL j (Ui + hUj^ + cUjL +• dU]^) + 

m-i , .I T - -ii r. . . . . . .i n "^yz z y 

1 d X 

Jadx 

- 

~ ® ”i- 

-\adx 

Then multiplying hy m^ e « and collecting terms, we obtain 

Ui -I- (a - ^ log mi \ Ui + hui +• cun +■ dui +. 

-Tcyz -^z -*xz X 

(28) 4 jc (a - ^logm^ )^ c^J u^^ 4- jh (a - ^ log m^ 

4^cL (a - ^logm^'^^g^ - m^J ^ ° * 

As before, designate the coefficients of (28) as follows; 



a - 



^ log mj^ 

^ X 



bi - b ^ 



Cl = 



c ) di* d ^ 



5l - c^a^ ^ b^a^^ + b^^ ^ » d^a^ ¥ ^ix " “l • 



Then (28) is in the same form as (l), and is identical with (25) 




•4 •$}■ » 
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Now, however, the x-invariants resiilting from the substitution 






~ ®1®1 “ ®1 = ®1®1 " ®-l®l ^ 1 ^ ~ 



a 

z 



- c 



A 

^ ^iog mj ^ 
^ z ^x 



ll^ - a^hi aiy - f 1 « ay - h^ - ^ 

“4 * ^1^1 + ^1^1^ + Cl^ly ■'- = 

= mj_ - dj. + ha + cay + ay^ - h log m3_ ^ 
- c ^^log m3_ _ 5"^log lOj^ 

^ 7^ 3C * 



If ]q^ » Ij^ • mjj^ "S 0, we may reduce as before, and solve. 
If, however, at least one of these invariants is not zero, we are again 
faced with the dilemma that our original hypothesis is not satisfied in 
the most general case, and the chain cannot he continued. Thus, we must 
again modify the hypothesis in order to iterate and continue the chain. 
Our hypotheses heccme 



(a) 1^1 ^ “ 0 ; 

(h) a s a(x) ^ h =r h(yj^z)^ c * c(y^^z)^ 
d = d(x) ^ e » e(x,y, z) — ac j f s-f(i,y, z) » ah ^ 



g * g(x) 




i 
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Under these hypotheses, m^jstad - g, is a fvmction of x only, as is 

^log How when we transform (l) into (25), the coefficients 

- - - - • 

3 X 

will he 

■ a - 3 log m^ , h2^ h , = c , d^ * d, 

— 

®1 * ^1®! f ^1 - ^l®-! > 8l “ <il®l •*"'3.12 " “-1- 

Thus we will find that 

'S 0 , l2|_ S ’O , mj^ = - d^ , 

hence, the hypotheses (a) and (h) are both satisfied once more, and 
the chain may he continued. 

Ohserve that, as before, as the chain progresses, the invariants 
must have the form 

^l+3i ^^l+3i~ “^1 + 3i * ^ =0,1,2,... . 

Thus if m]_ ^ 0, the invariants for the iteration will vanish only 

if m]_; is an Integer multiple of d^, i.e., ifor sane J»0,l,2,,, 

¥e also roniark again, that:.tp have a chain of y-invariants or s-invariants, 
hypotheses (a) and (h) must he modified accordingly, 

D. How we wish to consider methods for cascading the eq.uations when we 
employ the substitution (4), in the event that at least one of the 
^-invariants is not identically zero. For this chain we wish to solve 
(lO) for u in terms of u_p and, after the prescribed differentiations, 
substitute into (4) to obtain an equation of the form (25) in u_2^. 

In order to accomplish this we are again forced to place certain conditions 




I 



6 ? 



on the coefficients and the related invariants. Thus we introduce 



Method 3 ? For this we initially assume that 



(a) I2 " ^3 “ ° ♦ 

Solving (10) for u, we obtain 

u = — (u T + auT ) 



and hence 



““-V'' 

h 7 



«S » Jl (“-L + “-1^* ) - 1 + J 



-1 






y ^ z ^ -1 



) 



“^1 ^3^^ au.j^, (u.j^+ au.i^.+ ayU.^) 

-—2 (u-1 H- au_i). 

m-i aya^ ^ ^ ^ TT ^ 

Substituting these expressions into (!<■) and collecting terms we obtain 



(29) 4- + ^-1^-1:,2 ^ °“l"-lxy 



where a_i = a, b,= b - ^ log ^-l , c -i *»c - ^ log “*^-1 

••J- —X n . —X \ f 

d 7 d z 



p • 
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— d - c ^ log m_2_ - 1 ) ^ log m 3^ - ^^og m„i 

XT ” “ 



3 y 



+ ^ log m 2^ 

“T7— “XT^ ' 

^-l* ^-l2 + ®-l°-l f ^-1 *^-ly +^-lb-l 



and 



S-1^ a.iy^ +^-ia_lj^4 cj.a.2y+a_^d^^ 



- m 



-1 ^ 



axid (29) is of the desired form. Computing the new -x-in7arianta, we 
obtain 



and 



^ ■ ®y 



k - sr c — ' a ■ 
6 X z 



3^1og m_2^ 

a^^y ^ 

a\og “-1 .. 

“ 5 T 5 T “^ 



- m_2^ - 3 (c a log ^ h 3 log 3^108 

^ ^ 4 ' ' 






d ya 



- ^ ^08 ”^-l 3 log m.i ) - Sy (C - 3logm_i) 

a z 



% (t - 3 log m_3^) ^ 

d y 



As in paragraph C, we will he in a dilemma if the new - x- 
invariants are not all zero, fcrthen our initial hypothesis (a) will not 
he satisfied in general, hy (29), and our chain cannot continue. To 
escape this trap, we once more lir5)Ose additional restrictions on oxar coef- 
ficients, Our revised hypotheses will he 




i 
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(a) I2 ^3 =05 

(b) a = a(x), b -b(y,z), c ®c(y,z) , 

d =■ d{ 2 )^ e =-e(x,y,s) = a(x)c(y,z) , 
f =f(x,y^z) = a(x)b(y,z) , g =g(x) . 

Under hypotheses (a) and (b), m_]| becomes a fimction of x 

only, and hence the coefficients of (29) became 

s„l — ad - , 

With these coefficients, the -x~inveLriante for (29) 

^5 ^ ^6 ^ +• ( 2 ^ , 

which means our hypotheses (a) and (b) are again satisfied, and the 
chain may be continuedo As observed in methods 1, and 2, the invar- 
iant _ 2 j can vanish only if m„i is an integer multiple of 

that is, » ”J<ix» this point, that if m is a positive 

multiple of we should us© substitution (2), while if it is a riega- 
tive multiple of d^, substitution (4) would be the more advantageous. 

One more it is worthwhile to note that if we desire a chain of - y - 
invariants, or -z- invariants, hypotheses (a) and (b) must be modified 
accordingly. 

E. Finally we consider a method for cascading the equations when we employ 
the substitution (3), in the ©vent that at least one of the invariants 
in (8) is not zero, Althou^ this substitution appears the most natxiral 








rwipiiiii"" '• 4 ^^ ^ 






i. 




m 



v: 
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in that oar eq,uation (l) is iamediately broken down into three first 
order equations, we will find that the most stringent conditions on the 
coefficients are required in this case in order to generate the chain. 



Method 4; Oar initial assumption in this method is that 
(a) h2— kj^ ^^1 ^ must consider the system of equations 



U3. + 


au ^ 


"1, 


u^ 4 


bu^ = 


2 

’^l i 


2 , 


2 . 


2 


Ui 4 
■^z 


CUj — 


m^u , 



If we solve the second equation for u^, substitute this value 
into the first equation and solve for u, and then substitute the result- 
ing expression into the third equation, we obtain 



2 2 2 

Ui 4- cui = mie 

■^z 



-Caamf, Jadx-bdy r Jbdy g , 

[J® r 

;,z) j dx -h I(y,z)j ^ 



4 Y (x. 

where Y and X are arbitrary functionB of their respective arguments. 
Ve then proceed as before; 



Jadx r 2 sl 2 . 1 ■ 

m ~ ^ J ^ 



X. 



Taking the partial derivative of both sides with respect to x. 



( 30 ) 



^adx 



4 



^(a - ^log mi )(ui^ 4. cu^ ) 4 cu^ 4 



Jadx-bdy f Jbdy 

Me uj,dy 4 - Y 



=■ e 



j 



Wo then maltiplj- ‘botli sidos of 
derivatiTe with respect to jj to ©btaiai 
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ftd^-adx 

hy © ^ aad take the partial 



^ r(h- 3 logra|) r (a - (^1002^) (u| +cis|) + u^ -fcu^ 

“a? ^ >i y ^ ^ 3iz X J 



(31) 



+ (Oy - ^ 



Tr 

2\/,,2 __2 
'Z 



^ + 'x"V <^4 



^ \ 2 2 2 2 
+.eW]_) 4 (a - ^logm,3_)(u^ ®’-^lv'*"®y^l^ + 

b ^ yz 



0 






_2 



Fimlly, multiplying (31) hy laa^e'^ , and coUecttog terms, 
w© obtain an equation for Uj, which is of the same form as (l), hut 



2 

whose coefficients are = a - 



h| s* h - ^logn^i c^ — c ; 



■yr 



3T 



^1 * ^1®1 •*■ + ®L J ■*■ ®ly ' 

®1' + (4ei)y-^g?^ " m-1 . 

Computing the new invariants, w© find 

h^ = hi^ + hiei - di - - eiy ^ 

+ 4A - 4 ^ 4^ - 4xi 

>\= 4yi 44, - 

= 4 + “Vs ••■ (“M>s +4y4 +°i3, . 

We are faced with the same difficulty as before. In general, the 

invariants h^ and k|^ will not vardsh, hence it will b© impossible to 



2 
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continue the chain in this manner. We will he able to continue the chain 
if we revise our hypothesis to he (a) hg ^ ^ y 

(h) All coefficients are functions of 
y only, except c, and that c 
is a constant. 

With these hypotheses our new coefficients heccme 



^ = a; 1)2 1 - ^ . 

®1 = 44 -> ^ ®ly' 4 = - >4 ' 

new Invariants hecome 



I15 — \ - ll^ “ 

-}- csy . 

Thus the hypotheses (a) and (h) are satisfied once more, and the chain 

may he continued. We oh servo that in this method the chain will terminate 
2 

with all Invariants zero if and only if is a negative Integer multiple 

of cUy. 



These methods described in paragraplis C, D, and E are not, of 
course, the only methods by which chains of equations could he generated. 
We could do a method similar to method 1, say, with = 0 and =s m^^ 



or We could do a method similar to method 2 with kj_ the non-zero invariant^ 
or the non-zero invariant. There are many variations on the theme 

presented here, hut it is the theme itself which is the important thing. 

In any such method certain restrictions of the coefficients are necessary 
to iterate once, and other restrictions are necessary to iterate more than 



I 
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cmc©. The haBic idea is to find an e 3 cj)resBiaa for u in terms of the new 
dependent variable, which may or may not he an integral es^jressian. We 
then differentiate this ea 3 )ression as required hy the particvilar equation, 
i,e. the equation involving the invariant coefficients. Substituting the 
resulting expressions into this equation, we reduce it to an equation in 
the new depeodarct variable alone; if necessary, we separate these ex- 
pressions fraa the invariant coefficients and differentiate a^in suf- 
ficiently to remove any integral signs which may appear. If we can 
acccmpllsh this, we can in general Iterate the equation. 

F. We consider now the most general form of (l) vdilch can be reduced by 
these methods described in paragraphs C, D, and E, The question is, what 
must the coefficients of (l) be, in order that (l) may be reduced to a 
second order equation by these methods, when the x-invarlants or the . «*x- 
invariants are not all originally zero? We observe first, in methods 1, 2, 
and 3f that in order for the m-invarlants to vanish after n-1 iterations 
we must have 

g - ad = nd^ , where n =• ±1, ±2, ... 
or 

(32) d^ 4 ad - £ 



Equation (32) is Integrable, and integration gives us 



( 33 ) 



d s 




pl«a an arbitrary constant 



which we choose to be zero. 

Using (33) and the hypotheses of method 1, we see that our original 



equation must be 



7h 






r J|dx 



( 31 ^) 






-{■ g(x)(Wy+. 1».^+ U) s 0 , n^l^2, ,«,, 



for mathod 1 to bo applicable, ' Usiaag ( 33 ) aiid the fcypotbesos of mothoda 
2 and 3^ w soe that oar origiiial equatioca must bo 






/ -^fic f \ 

)^z +• b(y,z)t^2 + d(y,z)A^-t-^o \o g + 



( 35 ) 



^ a(x)c(y,z)Ti +a(x)b(y,z)W2 f s(x)ti so, n *1,2,,,, 



for methods 2 or 3 to bo applicable. 

Finally wo eraasider the fom of tho coofficioiats whoa method will 
bo applicable. In this case w© dbsorro that in order for tho m invariant 
to vanish after n-1 itoraticsrss, wo must have 



or 



( 36 ) 



g - abc 5SS. nea^ , n* 1,2, ,,, 



a-y + ^ 

n nc 



Equation (36) may bo intogratod to give 

C bdy ^ fbdy 
On \ o Jn 



( 37 ) 



1 



c n 



Using (37) and tho hypotheses of method 4 , wo see that our original oquatioa 
must be 



- Cbdy ^ (bdy . 

- Cbdy r Tbdy v 



+ cb(y)Uy + (c© 



c 
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+ 



, c n H 




for method U to he 



applicable , 



2 + s(7)u =0^ 



SKSTIOJS Y 



HIEERBOLIC EQUATIOKS OF N™ OHDER 
IK H IKDEPENDEHT YAEIABIES 



Ao Having considerod tlio extension of the Laplace cascade method to the 
equaticm of third order in three variables of the mixed derivatives type, 
we will now t\zm om" attention to a generalization ©f these results. That 
is, we consider an equation of n^^ order in n independent variables, 
in which only the mixed derivatives appear, and we shall attempt to show 
how these methods can be applied to such an equation. 

Before proceeding to the n^^ order equation we will first intro- 
duce an operator notation which will simplify soEowhat the invariant ex- 
pressiona and relations, which tend to become awkward as wo increase n. 

In order to do so, wo mention here a rule which has been devised by J, Lo 
19 

Roux and shall hereafter be referred to as Le Roux's rule; 

"For any differential operator D which is a polynomial in 
and the following relation holds; 






D(uv) Cilfl)(v) ^ ^ ^ ^ Dy(v)'^2] 



^x 



'dJ 



i 



D 



(v)^.-|^d;;^(v)1 













(v) t 



3^d'" (v) 



+ 3 



d"'(t) + 



^3u 

3y3 






3 



* • • 



where 





ft 

, etc. 
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This rule is an extension of Taylor’s formula, and is, of course, 
completely analogous to Leihnitz’ Rule for the n^^ partial derivative of 
a product u v. It can he similarly extended to products of more than two 
functions. As an example of this rule, consider the operator D(u) “ 

!U_ -h au^ 4- hvi„ + cu — 0, Here D = ) ^ o_a^ j.h^ 4 *c, hence 



D'=d + a , d's^xI) 5= I# "the identity operator. 

For this operator D, it is a simple matter to verify the identities 

(1) d’d'(u) = D(u) -I- hu 

y X , ; 

^2) DiDy(u) = D(u) ^ ku^ 

where h=Ss_ 4 -ah-c, k»h 4 -ah-c are the two Darhoux inviJpiants. 
^x y 

Then denoting u^=D^(u) and observing that D(u) —0, we have frm (l) 
(3) ®y^\^ * ♦ 

Operating on (3) with D^, we obtain 



(4) —Dj(liu). 

Using the Identity (2) to evaluate the left-hand side of (4) and Le Roux’s 
xnile to evaluate the ri^t-hand side, we find 




s: h Ut 4- ^ h u 

Jy 



=■ h Ut 4 ^h , h ,u 

^ 5 ? e 



= hu^ + 



D’(tt ) 

y ^ > 



or 




\ 



I 



\ 
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D(uq^) - ^logfa Dy (u^) + (k - h)uj^ = 0 , 

ay 

This may he vrltten in the fc*m "=■ 0 

where D-i — ^ j. a, ^ + h, 4- c, , 

d^U ^37 ^ 



aj^= a - 






glogh . 

*Ty 



= c - h ^logh ^h - a^, 

Ty 

and these coefficients agree exactly with those obtained hy Darhoux' 
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for the i^lrst equation of the cascade method. 

Our problem now is to expand these procedures to the n*^ order 
equation in n independent variables, of the mixed derivative type. Can 
we in this maimer predict the form of the "invariants”, the corresponding 
identities, and methods for cascading the equations? 



B. Let us, for convenience sake, change the notation sll^tly. Instead 

of indicating the independent variables by x, y, z, ...., we will find it 

th 

more advantageous to label them x^^, X2> • • . .2^0 Then let the n order 

t 

operator with which we will be concerned be defined by 



(5) 



D(u). 









where we always require that ij ^ ij^^ and that ijs only if 

ij^.lS 0 e The sum is to be taken over all acceptable eombinatians of the 
ij where lj=o,o<»o,ao In (5) we write Uj= a®*! define 

S Ue Finally we require that in all of those operators we will set 
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^12 ii ~ illTaetrat®, l»t n 2, Tliea 



D(u)sr + ^ol 4 ao2 4- 






For n» 3 w® 




A + * aju 

^3: ^3^ 



®^001 ^ ^ + ®002 + ®003 * 
^ *1 ^ Xg ^ 



and th® breTity of (5) "beccaaios mar© apparent with increasing n» Having 
seen these exa^les in explanation, there will he no loss of clarity if 
wo drop the zero s\ib scripts of the a’s wherever they may appear. Thus 
for n — 2, wo writ© D(u) = u^ 4- QgUg + -► 

In order to point the way, lot us briefly simmarizo the 2*^^ order 
results of farhoux and aur 3^^ order results in the li^t of this hofiatioK^ 
Zigr the second ordei* eq.uation we have the following operators; 






ail 

i„+a 

“2- 



+ “1 



a-L 'h 









The two Darhoux invariants are given by 




I 
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h 



= a 



^ a = D^(a^) - a ^ 

k r a ^2 + aj^^2' ^ 

whil 9 ths identities (l) and (2) tecom©, respectively 

D^j(u) * D(u)l liu = D(u) + D 2 (aj) - u ^ 

D;Ji) 2 (u) ~ D(u) 4- ku ^ D(u) + " ®' J 

For the 3^*^ order equation, the operators are as follows; 



^X3 



+ A- -I- ^2 iL + 

"1 



1 ^ r “■a -r- + *^3 4 — + 4 
PT 7^2 3^ 









.-2 

D,'U i 



12 



5^ 



■H ®'12 



= 4- + ^13 






_ ft ^ 

®23 - ^ + ®23 . 



and associated with this equation are the invariants 



y * 
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^3 == -^1 



2 If If , 

“l = ®X2“X3^“23>-^ 

4 - 



K = 



K =■ 



tt 



® 12 ^^ 23^“^2 

^23^®'12^“^2 



“2 - ^12^23^ ^13^*^ 

■4 = 



\ - ^13^^23^ *^3 

1 ■= d” (a ) -a 

2 23 13 3 



“1 = ^1 <^23^-^ 
“2 = 

“3 = 



It 



Dg (aj3>-a 



D3 (a.j^2^-a 



4 = 

— d” D* * (a ) -a 
3 23 13 12 



“-1" ^13(^1)-^ 

'• / V 

DT_(a^)-a 



-2 



= D^(a^)-a 



A typical identity involving these invariants is 



I II " '• ' 

Dj^D23(u) - D(u) 4. (D^(a22)-a2)u2+(Dl3(a23^“^3^^'*"^®l^®^23^“^^’^ 



C* Now let us generalize these results to the n^^ order linear differ- 
ential equation of the form (5)« In order to demonstrate the magnituiip of 
the problem, let us first compute H(n), the mmiber of different identities 
associated vlth the n^^ order equation, D(u) = 0, 

An identity for D(u) is a ccnibinatian of differentiated operators 
d!, d", such that successive application of these primed operators to a 
function u = u(x), where will yield D(u) plus a linear 

cooibination of the mixed derivatives of u of order less than n. We shall 
denote such a linear ccmbinaticca by H, for remainder, and define the coef- 
ficients of R to be the invariant coefficients (or simply invariants) 
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asaociated with the operator D. ¥e shall show presently that the order 
of R must he less than n-1, althou^ this is not evident a priori* 

To he a systensatic we proceed as follows: 

t 

Select arhitrarily a prlaed operator having order n-1, i.e. 
for some i«l,2, ,, «,n. There are, of com*se, n choices for this oper- 
ator. To obtain an identity as defined above, there is only one choice 

of operator such that +E. This operator must he 

(n— l) 

Dt o 4 n 4 I -1 siiwe this contains the differentiation which 

1 2««,i-l i+l«,.n ■ 

t 

is Hissing from D^. QJherefore we have exactly n identities of the form 

I (n-l) » . 

^1 ... i-1 i<fl...n^^J — D(u) + R, 



i*l,2,.. 



,n. 



Next select arhitrarily a primed operator of order n-2, D 

IJ 

seme 1 = 1,2, ...n and some J = l,2, ...,n, J?^l. There are n choices 

•• « 

for i, and n-l choices for J, hut since D total only 

ij <3^ 

n(nsi^} = (t) such operators, where denotes the appropriate binomial 
coefficient. If we wish to obtain an Identity when we apply such an oper- 
ator, we must apply it to one of three things; 

Sn-2) (u) 



1 ) 



'1 2 ... i-1 141... j-1 J+l...n 



2) (u) 

1 2 ...i-1 1+1. ..n 1 2 ... J-1 J+1 ...n' 

M (u) 

^1 2 ...J-1 J+1 ...n 1 2 ...i,l i+1 ...n . 



Therefore we have 3 ( 2 ) of the form. 





# ■ mm 



^ m 



1 
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L(u) s D(u) + E 

whore L(u) is on© of the throe oxprossion shovm aboro. 

If wo continue to count the identities in this manner^ we find 

that 



( 6 ) 

where 

( 7 ) 



n-l 

H(n) • XT “j (j) > 



01 . 



- 5 :) . ■ 

d 

. ’ IvV sot of all ordered, integral partitions of J, 

To verify this formula it is only nocosary to count, as wo have 
already done for the terms J— 1 and J = 2, the identities for the J 
term. For each of these identities, the last operator to ho applied is 
. , where the i fcsrm all the possible combinations of n 

IjXg o • • Xj p 

integers taken J at a time. Hence there are (^) choice for this oper- 
ator. To comploto an identity wo must apply this operator to a combina- 
tion of operators "vdiich contains each of the differentiations . ^ , 

P 

p = 1,2, J, oxico and only once. Thus, for example, w© may have a 

combination of J operators, each of which has one differentiation ^ 

or wo may have a combination of J-1 operators, on* of which contains 

^ 2 

the settcnd order differeostiation O while the others eonta-ini’ 

3 iVik. 

ip H 

only first order differentiations ; and so forth* The miaiber of 



cholcoB for oach operator can only to tlie number of ccanbinations of the 
differentiations still to be included, taken as many at a time as the 
order of the operator being considered. Finally, since the order in 
vhlch these operators are displayed is a factor in determining which 
identity we have, wo must consider an ordered integral partitions of J . 
The summation in (6) is taken only up to n-1 since for J —n, we have 
the trivial identity D(u)=D(u), and this is not to be counted in our 
dot eminat ions. If wo tabulate the appropriate values for n=2,3,it-, and 
5, wo have 



D. We turn now to a consideration of the invariant character of the 
coefficients of an arbitrary R. Consider first the identities of the 
form 




R(2) =■ 2 




N(3) = 12 




(n-1) 

®i 2 ...i-1 i-fi... 



^ (u) — D(u) 4- R. 




• • J-1 J+1- ■ • » Sj-^.iKi-*'4-i^J«' •• 



Cji-3) 



j^l 

n n 

|i fci - ,j.^ 

j^i X^L ^ ^ ‘ 



-f 



i^l' ‘ ‘^- 3 ?V]‘^ 



+ 



We may write tMs more ccmpactly as 



V 



(8) B^(u)iuj“)-^V 



n 



J^l 



. .(n-2) . 1 
^J) (IJ) 2 



n 



• • 

jjti 



K*l 






+ a^\f. 



j/K. 

To Ulimiinate the notatioa^ the superscript on u indicates the order of 
differentiation, the subscript on u, which is in parentheses, indicates 
the variables with respect to \diich u is not differentiated and the 
subscripts of the a coefficients, which are in parentheses. Indicate 
the integer or integers emitted from, the subscripts of the a coefficients 
in D. For example. 



~ ^12 ... J-1 J+1...H 



U) 

^(jk)^^12 J • o k**l k4*X*o«n« 

Observe that factor for similar reasons, the 

a term in this expression contains the factor 1 . 

If we scrutinize equation (8) -we find that there is ^ust one term 
of order n-1, and its coefficient is 1; there are n-1 terms of order 
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n-2, ajad th«ir coeffieieats ar® th© (^2^) 't®™s of 

orcL»r (a-3);» with eo«ffi©i®nts aaad .in g®n®ral th®r® ar® 

t®rms of ord®r n*»q.* Consider now th® ®ff®ct of th® operator 

^ 

^l...i-l i+1 ... a ^ '^'‘(i)- 



-h a-Qy, ■vh®n it is applied to D^(n),. 



J#i a *i 



(9) 



(jffi-2) I 'v 

^ ^ ®'(jk)^(jk) *** >r^ 

Li , 

J¥k 



Equation (9) tells us immediately that there are no in'variant coefficients 
for the terms of order n-1. This statement is readily seen to he true 
no matter '«hich identity we consider. Thus a discussion of the invariants 
must begin with the terms of order n=2, 

Fran the fact that the indices i, and k are distinct we con- 
clude that further note that (9) contains 

(n-l) -h (n-l)(n-2) different mixed derivatives of order n-2, and */ 

these are all th# possible mixed derivatives of this order. Thus in order 
to obtain the identity 



1+1 ... n “!(“)=»(“) + 

l.e. to make the expression D(u) appear on the ri^t hand side of (9)^ we 

must (at the very least) add and suhstract to the ri^t hand side of (9) 

(a-2) 

*(«)“(«)• 



n«2: 



This results in the Invariant coefficient of order 
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( 10 ) 



.t ^(j) •ur’ “(ji) ’ 

Sh 



Sixico there are n choices for i, ajod n-1 choices for we have 
a total of n((n-l) invariant coefficients for the terns of order n-2. 
We will now prove the following; 

Theorem 7 ; The invariant coeffielenfas for the terms of order 

n-2 which appear in the reiminiler R of any Identity are of the form 

5^ . 

These ”lnvarlants^are all true invariants, l»e. invariant tmder the 
change of variahles u — X u, and no Invariant coefficient for terms of 

M U 

order less than n-2 can^ln general^he a .true invariant. 

Consider any identity given hy a sequence of operators. 



(a) 

e o o - 






Dl ' , Pt ' 1 (u) =D(u) R. 



; p(=) 

,ig^ ^***^q^ ^1***\ 

where (n-a ) ■+■ (n-h) + . « . -J- (n-s) — n. 

If we write this out in more detail, we find that such an identity is 



actually 






(n-a- a) 



— CL. 1 -{■ d + »< 



Yn-b) 



^(n-b-i) 






ct ^ ^ ^ *f* 

(i^ Jl 5 ' ' 'jJ*b) ^ Jr ' * Jb> 

\(tI”S) n v(n-s-l) 

^ ^ , + y a . d ^ u 4 . 

^Xj.* • • <^Xj^ (j) 

i^s) j 



(1^1 



• # # 



= D(u) R . 





p 








I 



I 




I 




A 



■■SI, 



* 
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By aa jLnsp»cticm of th® l«ft hand eid® of (ll) w© can determln® how ©ach 
t©m in D(h) is aiad© to appear on th® right hand aid® of (ll) . To 
obtain th© n^ ord*r term g w® laost ccnibin® th# hig^st ord®r 

, ani apply 

il**«ia, «3l»*«^, ^l**®^q 

(rh® dotormination of this n^^ order 



differentiation in each operator D 



th® rosnlt to U 

o o V©"'3 
(i-:i l]_«««;>lg) 

term is, in fact, th® criteria used to detenain® -which ecmbination of 
operators yields an identity*) Th® terms of order n-1 which arise are 
obtained in two -wayso They are obtained by combining -fch® highest order 
differentiation in each operator Di*l..ig^, appiy= 

ing the result to each term of the sum 

n 



(12) 



U) 



(ss. « S-l) 



r 

<3 1 * 3 x^» ® * 

(i,j »-l^,..lg5 j) 



Using ^Leibnitz's riile we see that, among o"thers, the terms 



“(j) “‘ut’ 



J “ o O 



must appear. Such terms of order (n-l) also arise when we cooiblno, one 

by one, a differentiation of second hipest order in any one operator 

i n Tr with the highest order differentiation in erery other 

^l**°^a ***^ ^1*** a 

operator i j*s* I>l^»ook s applj this result to 

® a • ^ 
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(13) 






Tlio coefficients of the terms of order n-2, ■R&ich are the terms 
ve are coaacemod with in this theorem, arise in four ways. We may ob- 
tain siich terms by ccmbining the highest order differentiation in each 
operator applying this restilt to each term 

of the sum (12). Again using Leibnitz's rule, we see that, among others, 
the terms 



9 ^(4) . u i,J = 1,2,... n; lit J; 

i,J^l^...l^^ 

must appear. If we combine, one by one a differentiation of second 
hipest order in any one operator . , with the hipest 

K2^...Kq^ 

order differentiation in oTery other operator D^^) . - *1.1©^*?:) ’ and 

apply this result to each term of the sum. (12), we see that, among others, 
the terms 

(n— 2) 

i, J —1, 2, ...,nj i ^ J 

1.J 



(n-2) 

^(j) ^i) ^(iJ) 



will arise. In a similar manner, if we ecmbine a differentiation of third 

hipest order in any one operator with the hipest 

order differentiation in every other operator ^PPly 

this result to (13), w© obtain the terms of the form 

(n-2) 

®(JJ=) ""(Jk) ' 

while if we ccmbin© the hipest order differentiation of each operator 





*■ 
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V t Q'PPly tMs result to oac2i term of tlie sm 

a (n-s-2) 

^ „ U . 

(l-,J:,k •1^,.,1 j 
s 



we must again obtain terms of the form 

(n-2) 

^Jk) "(Jk) . 



Comparing the coefficients of these terms of order (n-2) and recalling 

' e 

the procedure necessary to make D(u) appear on the right hand side of 

(9) > we see immediately that the invariant coefficients for the terms 

which appear in the remainder S must be of the prescribed form 

( 10 ) . 



Knowing that these Invariants must be of this form, we proceed to 
show that they are tni® invariants. We consider the equation D(u) = 0 
and the change of variables u - Xu > where X(x) is not zero. Using 
the chain rule for differentiation of a product, and employing the notation 
of (8) we have 



(14) 



n 



n n 



(n) 



u 



-f 



Z ^x 



-(n-1) 
"(i) ^ 



1 X X ^^x 

2 i*l j*l Xj 



- (n-2) 

^(ij) *'’••• 



\-(n-l) ^ \v «(n^) ^ ^ ^2 X „(n-3) j. 

•A^ (i) + ^ ^ i Z 2_ Jr., A, u 

J=1 T^. 2 J3. k»l <5z.(3x^ (ijk) 

J^i ^ jel k#l * 

i«l,2,...,n 

J 



lu) 6 Ts: (ww 



1.2. 



kij 

iiH 
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and Bimllar other expreseions involTiag derivatives of order leas than 
(n-2)o Substituting the relations into the equation D(u) » 0, and 

dividing throu^ "^7 ^ f we obtain 



5(S) = 




n 



I 



^ logX ) u 

Tx| 



(n—l) 

U) 






n n 



Let us indicate these nev coefficients with the obvious notation 



^(i) - H±) 



^ logX 



“(«) + *(!) 

Then the new invariant coefficients (associated with D), for the teims 



u, . , are 

(ij)' 



and In teinas of the original coefficients, these become 



“ a. 





^(J) 




as was asserted. 



■ 

I 

I 




92 



To show that ao laTariamt sooffiei«at for terms of lesser order is a true 
ia'mriaat, we must first laTostigate further the character of these other 
iamriaats* 



E* ¥e will now giT'o a general method, for the deteminatioa of the in** 
■Tarlant coefficient of any order term in an arbitrary remainder B, 
Consider any identity given by a sequence of operators. 

V**^a ^i*“ q !••• s 

where (n-a) +• (n-b)-f ••• +('n-s) = n. 

Wo will denote for brevity 



S^) 



... D. 



(a) 



“J- <l 



-<0. 



Suppose wo wish to find the invariant coefficient for the derivative 

„(p) when c is a term In the expression <i (ii)« 

m....m^ ' o mj^...m il**’-le 



S^nce each differentiation ^ , i = 1, . • .n, never appears in more than 

^x. 

one operator on the left hand side of (15)> wo see that the left hand 



side of (15) contains the term 

0(Co^ nr > 

and no other term on the left hand side of (15) contains u 



(P) 



HL-i • • #in. 

1 p 

a factor. Thtis, recalling the procedure necessary to make D(u) appear 

on the ri^t hand side of (15)# we see that the Invariant coefficient 
(p) 

for u must be 

m^...m^ 



as 



"O'!'.) 



(16) 












4 
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Wow suppose we wish, to fiiid tJae invariant coefficient for the 



.(p) 



derivative 
which contains u 



when there is no term in the expression (u) 

^ as a factor. We may assume that i (u), 

... .. ■ (p-j) 

contains a term Cq ^ for sane ...p, \diere we set 

(o) 

u. This statement is obviously true for, at worst, the expression 
(s) 

•••] must contain a term c^u, (although ®ay vanish identically). 

We may then apply LeRoux’s rule, extended to an operator with n inde- 
pendent variables, to the product c . Since the original 

operator D never contains any differentiation a multiply, for any 

Jx^ 

1 »1,2, n, the same is true of the operator ^0" . Thtis our extension of 
LeRoui’s rule for this type of operator is 



(17) B(uv)=»d(t)+|;_ It •••' 

J-1 e* J,k=l 

and hence 



. , jy(l> , , <P) 4. 

®»p-JU!::®p ® mii».mp 

From this wo easily conclude that the invariant coefficient for 
must he 






( 19 ) ^ 



U) 



Vj+i'-“p 



(cq) - a 






a o oHl 



I 



9h 



Wo nsuLst point oiat tliat fomula (19) does not giro tho camplot© 

pietvsro. For if ^0 is tho idoatity oporator, then it must Tt>© 

VJ+1**‘“p /X 

that that . .nip appears on the left hand side 

of (15) • That is to say, there is no need to add and suhstract 
(p) 

ri^t hand side of (15)* hone© there is no term 

containing in R, In other words, there will he no inTariant 

coefficient for u^^^ • The same conclusion holds whenever 

(a) 

a appears as a coefficient in any on© of the operators D 

m^.»»nip 



ii-V 



x0>) 



x(<l) 



. This is apparent from tho loft hand side equation 



(D), and tho discussion which follows it regarding the determination of each 

term in D(u) which appears on tho ri^t hand side of (!!)• To summarize, 

whenever a _ appears as a coefficient in any operator , 

(h) (s) , » 

D. 1, which is on the left hand side of (15)» there 

Jl****3h - 

will he no invariant coefficient for u'^' appearing in R on the 

V-^p 

ri^t hand side of (15) • Otherwise, there will he an invariant coefficient 
(which may he zero) appearing in R on the ri^t hand side of (15)# which 
will he given hy formula (16) or formula (19) as appropriate. 



F, As yet w© have said nothing about the character or number of quasi - 
invariants, i.e. the Invariant coefficients for derivatives of u of order 
less than n-^» As we -have pointed out, however, there will ho no invariant 
coefficient for a derivative u^®^ when the ^ coefficient a appears 

m2^..»iap 

on the loft hand side of any such Identity (15) • From this fact wo may 
easily compute the maximum nunsher of quasi -invariants associated with any 
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identity (15). 



3Si» total mmiber of terms in D(u) wMcli are of order less tlian 



n-2, provided no coefficient is zero, is simply ^ fi 1 * sub- 

tract frcoi this the total number of coefficients of derivatives of order 
lese than -n-2 in D(u), which appear on the left hand side of (15) > 
assuming all are not zero, we will-have the maYiTmiTn number of quasi- 
invar iantS;>which c«n appear in E, on the ri^t hand side of (15) « For 
if such a coefficient would appear on the loft hand side of (15)# were 
it not zero, then, regardless of the fact that it is zero, there will he 
no quasi-invariant cctrre spending to this term. This number, then, is 



zero. For n»3# we observe that an unique qtiasl -invariant appears as the 
coefficient of u in each identity, hence there are a total of 12 quasi- 
invariants for the entire third order case. 

For n>3, however, the picture beccanes somewhat clouded. This 
occurs because while no quasi -Invariant can appear twice in the same 
identity, the same quasi-invariant can and does appear in more than one 
identity. Thus the determination of the total number of quasi -invariants 
for a given order n seems to be a rather tedioxis, lengthy, and unre- 
warding problem. To obtain an indication of the enormity of the task in- 
volved, however, it is quite simple to obtain upper and lower bounds on 
this number, which we shall designate M(n), 



( 20 ) 



(n-a)f (i 

For n 2, the nusiber of quasi-invariants is obviotfeiy equal to 





I . 
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An Immediat© lower ‘boimd is a consdgxsaace of tli® fact that an 
unique quaei-inTariaat results in. each identity as the coefficient of Uo 
Hence H(m) is a lower hound for M(a) o To obtain an upper hound, wo 
merely ohsoi^e that the total number of quasi "invariants associated with 
D must he loss than the nuniber which woxald exist if every quasi"invariaBt 
appeared in only one ideatity® Now with every ordered integral partition 
of Jji r I ^ there are associated 




identities of th® form’ 



(21) D(u) 4 - E, (see 



eqx»tions(6) and (7 ))o 

With each of these identities are associated a maximum number of quasi" 
invariants, given hy (20), namely 




Thus if wo auitt over (J\ , the sot of all ordered integral part it ions 

d 

of J, and then sxm over J, we will obtain the nunibor of quaBi"invariants 
which would exist if ©very qtiasi -invariant appeared in only one idesfcity* 
This number is 



( 23 ) 




where 




I 

I 



97 



Then for n ^ 3, we have 




N(n) ^ M(n) ^ H*(n)« 

Ccmputlng a few values for M* (n) w© find the following: 

n K(n) M* (n) 

3 12 12 

h 7^ 362 

5 5^0 '8510 . 



0 . Now we are in a position to coaig)lete the proof of Thooroni V. We must 
show that no invariant coefficient for terms of order less than n-2 can 
in general he a true invariant. We shall prove this rigorously only for 
identities in which/^ — ^ 4. and shall only indicate the proof 



for other identities. 






Consider the invariant coefficient for u , where p<m- 2 , 

p mi,,,mp 

and m. r for any i 1 . 2 , ,,,,p. Then c.u is contained 

1 ^ 1.0 ^ mi,,.mp.i 

in D*(u) and in fact 
r 



(25) 



c ■* a „ -r 

o m^,«, ^ ^ j 

Cn *. a ™ r 

^ “l*** “s>-i^^ 



where of course r must ho properly ordered with m^, ,,,,mp 







r 
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How caasider th.9 ehaaga of Tarialiles u “ X u, where A-X^^^ ^ Oc 
We extend the eqmticrag (lit), to obtain 






" .m + — 



N n-1 
o n 

(26) X ^2^* • ^ 



^n-p-1 ^ 

... ^... <»=^' V-**“p 

pa pa-i p+i-a n 



• • « 



o » i ^n-p 



= \« +... 

^ “i-"s 



If wo substitute (26) into D(u) » 0, and divide throu^ by X j. 
we obtain D^(u) » The coefficients in which wo are interested are then 



n 






.lOj, = . .Kj, + • - Vj ^ 



Jap-H 



(27) 



n 



1 + 



j*K. n 

+ Zi 



I i» a o < 






j,i:,i-p*i "1— ^jv*! x^vTs^ 



... 1 
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a 






n 



J.pfl 



, . .m m r 
P J 



1 jx 
X7^ 



n 



-^- V + . . . 4- 

j,krp.i ^ 

/V J “k 

J#k 
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A o o Yl-- “X 



whore 



\ 



=J^ ; 



"“(r) ' ®(r) + 4 i^- 

A 

In all tho ahove we must ccimmont that the suheeripts involved must all 
he properly ordered In accordance with tho conditions of (5). 

The quasi-invariant for Uja ,,,ia this particular identity is 

!*•• p 

simply Q = /T(c^) - , or 

Qa ^ (a „ V.) + « «) - ^n, 

(rj m^oeoDlpr inQ^^eein^ ^ 



Ihen the transformed quasi -invariant, for u , vill h© 



(28) Q S ^ (a 



) 4 ” (a 

*g- ^ HI • • •m r ' ziu • ♦ • 

^ IP (r) “i p 



) - a 

. .«mr' m^, ..m 



P • 



Substituting relations (27) , we find, after a sanewhat tedious computation 



that 
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n 



A ^ Illo 



»)^p+^ 



n 2 ^ 

■2L(^ N—VjV’'^'‘W°“i— “p“A’' ■^— Vj’O T3^x ^ 

l.lr=i*nal ^ /a 



J,k»p4l 

J^k 

J^k#p4^ 



(29) 



+ . . . 4 " 



n 



■21 (|r ‘“(J) 



a/.v) + ®-/„\a 



(r)“(j) 






Wo obsorro that th® coefficient of each derivative of X in this 
expression is itself a quasi -invariant vhich arises frcm the particular 
identity we have chosen. 

Then using an obvious notation we write (29) as 



Q »Q 4- T Q 1 iX_ 

(30) 



H- z: 1 



A5^ 



•I" 




Now if p*n-2, then Q 



^ ^ i 






o 








I 
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But 



'‘n=4 






.n 



1. . .(r-1) (r+l) . . .n 






(1) + .1 - = 0, h»«:e 

^ r 



Q * Q as was proTod previously in the moro general case. 

If p>n-2, however, then "q = Q -4-a linear combination of other invar- 
iants arising fran the same identity. Thus IT ^ Q iJi general, for the 

T t 

particular identities D_ ^ _ D (u)»D(u) + R, 

i r-i r+i ^ 

Althou^ wo shall not prove it rigorously here, it easy to see that 
for p>n-2, no matter which identity wo choose, terms involving 

^2 \ 

,—2 ! — d , .,,, will arise in the oa^ression fcr any transformed q.uasi- 

invariant, since these factors appear in the transformed equation D(u)»0, 
Thus, wo assort that for any q.uasi-invariant Q, Q Q, and in fact wo may 
expect that Q*Q + a linear conbination of other invariants arising from 
the same identity. This is the justification for the name "quasi -invariant, 
and this ccacpletes the proof of Theorem V, 



H, In order to see the problems involved in applying the cascade method 
for the n^^ order operators we have been considering, lot us consider a 
particular kind of identity once again. Wo will study identities of the 
form 

(31) ^i*»»r-l r+lo*»n^^5 " 
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n 



[ 31 ) 







whore Q 



li • • • 1 



1 # « « ^ ^ 
n-3 

ipfr, p=l,...n-3, 



is the auasi-iavariaat associatod with .u 






Donoto , (u) = ^u^a/ \u = u' 

1...3>lr+1...3ii YV “ ^ > 



3^ 






~ ^^(r) ^®'(r) ^(j) ' ** Jr) obsorriasg that D(u)s0, wo writ® 



(32) 



I 

u 



3“.+ “(r) (“> 




^J^(jr) 



i" 4 ^4 1 

^•••^n-3 ^1*** n-3 



It would bo most desirable if all the T's ajod all the Q's in (32) were 
identically zero, far if this were true we would have succeeded in reducing 
the order of the original eq.uation by one, resulting in a single first 
order equation, coupled with an equation of order (n-l). Ideally by this 
method we would like to reduce the n^^ order eq\iatlon to n first order 
equations, but this happy result will, unfortunately, occur very seldom. 

In the more probable oTont that not all the T’s and Q's yanish 
Identically, we would then hope at least to be able to cascade the equa- 
tions in the manner of the preceding chapter^,: mtU such time as the T's 
and the Q's may all bo identically zero. To this end it will obviously 
be necessary to lrq)oso certain rather rigid conditions upon the coefficients 
of D(u), Lot the hypotheses be 
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(a) T^= '“^- 3 ”'^ ^ ° ' 

(b) all coefficients are fvmctiona of only# 

(c) *11 coefficients with an ”r** subscript are identical. 



Then 








(33) 




')=T ^ U4 i 

^ n-2 


, l/=r. 


Solving 


f 

U s 


3 u -j- a, .u , for u. 


we find 



use 
hence 
u 



Je " 4=V ^ *'1J 



^1* * *^n-2 



= > L 









- e 



tW^r ^ , r 


f >5v)^ u'ta + F • 


^=^1 ...3^1 *' 


J J 



where q is the nuiaber of non-zero i^. 

Siibstituting these expressions into (33) we obtain 

1 n-2 

We then take the partial derivative of both sides of this expression with 
respect to x^, to obtain 



\a/ vdx 

0 ( 3 ^) r 



e 



[a(r)I^(u’) (d;(u*) . ^lorfr lC(u )1- 

L ^ "3”^ -* 

= 3i„. 






ip^r . 






Cancelling out the exponentials and collecting terms we find that we have 
an expression 

(u>) :r 0 , 

where is an n order linear operator of the form (5)^ but with 
coefficients differing from D. These coefficients we may compute to be 



^'(r) = ®(r)- 



logT 



/V 



^1 In 1 ^n 



where i =r for some p-1, .,.n, 
P 



a 



i. ...i - %.c.i ^ ^ 

^ n-2 - n-2 + 



\ ^ ^ ♦ a' “ T , where 

1... n-2 (r) 



^ X 



/V 



ip^r for any p=l 5 .«,., n-2 <> 

If now we compute the new invariants we find 




*(j) 

3 



j =1,2, ...r-ls r+1, ...n 



and Q' , 

^1 • * 



= T - 




rv 



But by hypothesis (c), a/’^'\ = Ec „ r for all r and all perrauta- 

tions ii»««in -3 » ^p^ ^ n-3* 

Thus if the new invariants are all zero , we may reduce the order of the 
equation ^D(u’ ) = 0 by one. If the new invariants are not all zero, the 
hypotheses of this method are again satisfied, and we may cascade the 
equations in the hope that the invariants will eventually become all zero. 

This can only hai^pen if T = n ^ f ^ is some positive integer, 

^ Xj, 
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other methods can be doreloped for identities of this form and for 
identities of other forms, but in every case, severe restrictions must "be 
placed on the coefficients before cascading may begin. We see then that 
the problem of attaelcing these eq.uations in this manner becomes extremely 
difficult and restrictive as n increases. Far n>4 the problem becomes 
enormous, with staggering numbers of invariants and ideitfcltios to be manipul- 
ated, and a large number of restrictions on the coefficients necessary to 

cascade the eg.uatlons* If it is desired to cascade the eq.uatlons it is 

1 

recommended that only first order, substitutions be employed. 



1 











I 



I 



I 



,^m 




SECTION 71 



A BRIEF SUMMARY OF EARLIER EXTENSIONS 



A. Extcnsloa of the Laplace cascade method have been made hy maisy not« 

able mathematiciams. One of the earliest of such extensions ¥as laade 

21 

by Darhoux himself. Darhoux considers the following system of equa- 
tions of second order, with n independent variables. 

Lot ^ system of n independent variables * 

and consider the system of equations , 

(1) ^ ^ i i,ka0,l,2,...n-lj i#k , 

This system contains n(n-l) linear hyperbolic equations of second order. 

2 

and we observe that this is an ovordetermined system, whereas the system 
considered by us in Section III is exactly determined. Darboux seeks 
first a necessary condition to give n linearly iadepondont solutions 
in addition to the trivial solution u« constant. To acccmpllsh this, 
Darboux forms the third derivative by taking ^ of both sides of (l). 

Interchanging the indices 1 and k, and equating coefficients of like 
derivatives of u, we obtain the relations 

(2) - ®'l'^k ^ ^-ll^ik ■ ®ik^lk (l«fe^i). 

Interchanging the indices i and i does not change the rl^t hand side 
of ( 2 ) and hence 



"ik 5 






Ik 
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Holding k constant, wo so© from Intogratillty conditions tliat thoro 
must ©list a function, call it lo^, such that 

(3) " -f- ^ , (11*1!:), 

Than condition (2) hocomos 

(It) I ^ _i 

and th© systom (l) has the form 

(5) ^ "5“ - _i. 3 .5 U-. 0 , 

% d?i 3^^ =1 Sft ^ 

k ■“ 0^ 1^ 2 j • o Ji**X j i^ko 

-J. 5Eic^ . l_ 5Hi 3u 
^ 3plc ®1 "5^1 

Usixig condition ( 4 ) it is ©asy to Torify that 

(6) f^j^(H^)=0 ^ ii^k^l. 

To start th© Laplac© mothod w© assume that v© hay© an ©quatlon 
of th© form ( 5 ), satisfying th© conditions of intograhility ( 4 )« Thus 
w© considor th© particular ©quation (for fixod i, k), 

(T) fii(u) - 0 . 

Darhoux than considers th© substitution function W defined by th© 



relation 





"‘’7 
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( 8 ) 



^ = -L- 



If VO substitute (8) into (7)^ tbo equation then beccmes 



i^yk ./dio^v . 1 H. \ , 

( 5 ) <-''>=<> 



¥o see that our Invariant h_., is nov 

ik 



^ik® 



^ - 1 3=1 J, iScV 
^ H], a^j-’ 




and if h^^" 0, the oq.uatian (9) has the sln 5 )lo solution 

e constant ^ 

We my gain a little insight into these proceedings if ve adopt 
the notation of Section V. Thus, for oaxiation (7), 



ik 



and 



ik 



Pi 





- 1 ^ - 1 b — 


k 


^i ^ ^1 ^ 






“ ^i 





Then the substitution function ^ is defined by 



( 10 ) 






vhilo the invariant h^^ W to written 
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(11) 




hr = 

ik 


D (a ), 

ik0^ ik ^ 


and hence 


(9) 


becomes 




(12) 


/ 

^k 

a 


(••V. 


(u)^ - ^ik^^^ "^(^ik^ ("®ik)^^ , 



If D (-a )=0, wo BOO that (10) reduces the order of 

ik ' 

( 12 ) hy one, aixd the equation can he solved by quadratures. 

In the more likely event that ^ik^ ^ cascade the 

equations in the following manned • Let 



and 



Then 









and satisfies the equation 

( 13 ) - 1. ^ ^ i ^ 0 

' ■ ■ 'p ■ 

In like manner, if we introduce the quantities Lj^» . defined by the 



relations 



no 



thou the fimetio* \J , for distinct values of i* and k’, satisfies 
the system of equations 

(IM = 1 + 1 

V ^ V 3^ ^ 

vhich is of the saioe fom as (5) • -» Hence we may iterate the suhstituion, 
and cascade the equations until such time as the corre spending invariants 
may vanish. 

( 22 ) 

B, In 1899, J. Le Roux extended the method of Laplace to linear 
partial differential equations of greater order than second. Le Roux, 
however considered the oqwtion with only two independent varlahles, 
while we, in Section 7, have considered the equation in n independent 
variahles. In the notation of Le Roiix, the equation to he considered is 
the order equation 






(15) D(z) = 

where n, n, n. 

If we suppose that the highest order of different iat ion with respect to 
X in equation (15) is equal to n- p, then the collection of terms 
which contain such a differentiation are 



Nn-p^i n-p+i-1 

^ ^ z L -jj ^ z -f* ••• S 

- n-PJ 






a-p 






T 







I 



Ill 



Wo denote the oxprossloa 



Jyi 57^1 



**“S 



^n-p 



as the differential multiplier of the term v ^ « and wo further 
assume that the coefficient a is equal to one* Le Hoiix considers the 
equation (15) from the point of view that there exists a particular 
integral of the form of Euler 



(16) 






V- 



where X is an "arhltrary' function of x, and the coefficients u 
^e functions of x and 

If we regard D(z) as a polynomial in ^ and , we may introduce 
the notation of Section V. Wo designate /t> 21s expression oh- 

tained when w© differentiate D with respect to .3 , p times, with 
respect to » q. times, and apply this operator to the variable Z. For 



instance, the (p^^q) 



th 



derivative of the term 




is 






oc(o<-l) •• ■ (oc - p n) g (g -l) " -C^- n tf) j y — J 



and hence 



(17) D^^'^^(z)in(n-l)...(R-P- 4 +l) 
^ y 



K'PV0-<1 







I 
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Since (15) la aBsusasd to have order (n=p) with respect to i, we see 
that the coefficients of (l6) mist satisfy the eq.uations 

1 _(a»p) 



(18) 



1 

(n-p)j 



(ir55T w ® 

t ^ j. (n-p-1) 

(n-p-1) I ^ 



(n-p) 



(n-p-1) 



^ ■^Tiiorr v-p-i + 



In the case when the variable x is a simple characteristic variable, Le 

Eoux considers the transformation 

(n-1) 



1 “ ^ D^-1 ( z) 

ins ar 



^ z 4- nA , z, 

3T 



and defines the functions and , 



— z — ze 
" ^o ” 

(19) 2 Hn-1^ , U - 

’ l':57. 



We then consider the eq.uation resulting froa (15) by means of the trans- 
formation (l9)j a^d collecting on the ri^t all terms which do not contain 
any differentiation with respect to y, we obtain 

(20) A(w'i')=X<, |^+ 

where A (t>^. ) designates a differential expressim of order (n-l), in 
which the coefficient of ^ ^/px is equal to unity. The coefficients 
are given by the formula, 

(21) 






P.P-l(Uo) = ^ ^P.P-l(«o). 
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The order p of the rl^t head side of (20) is, in general, eq.ual to 
(n-2). It will he less only if 

(n-2), . 

D^-2 ° • 

Le Eoux then proceeds to define the invariants h, hy the follow- 



ing, mile. If the 
independent, then 



( 22 ) 






X’s regarded as functions of y ar0 all linearly 

Xo ^1 



>^i 

• o o J. 



3 y 

• • • 



^ Xx • • • ^Xi 
I? JT 

o e 3 e o e 

}\i ... 



which are (pfl) in nmiber, setting If the X’s, regarded as 

functions of y, are not linearly independent, we define X^f^ as the ' 
first A which is lineardy Independent of \q* &a the 'f'll?et 



ils linear^ 


laidepen^ent' < 


Df Xp 


is replaced 


by 


Xo 


X«fi ■ 








SAo 


^XoCj, 




(22 ») ^i 


= 


T7 


yr 


JT 






e • • 

^Xo 


e. • 


0 • • 

^ Xe<i 








b y^” 


ay" 



O' ' '2 

c. Then the detemlnant (22) 



To show that these are generalizations of the Darbour invariants, we must 



consider the transfonnation 



.2 * 




z' f(x,y)^ 

<?(x‘) 

^ (x*,y*) . 



( 23 ) 



The first of these transforiGaticiris is effected without changing the value 
of , siaise u' s , and hence 



u' 



u. 



The transfcormtion of the indepeaadent variables, however, causes eq.uation 
(20) to heccsae 

(24) K'oiiL*K 



where 



{J* - U ^ ’ 
1 1 iP 



^ a ;* 



n-l-T) 



If- 



(25) 



4 . 37 4 10 

^y^ / 



^^ 2^®21 ^ i "^20 ^ o ) 



U ^ 1 

#00 



the coefficients © being dependent only on x. If we use the relations 
(25) to campute the new invariants h^ ; w© will find that 






\(i+l)(i+2) ^ (n-2-p)(i+l) 

2 <7 X 

“3x^ 



Since in general p»n-2, the above relation, reduces to 






(i^l)(i»2)_ 



h 



1 ' 



Thus the detsrmlisant h. is reproduced mutllplied by the (i4-l) 

2 
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power of the Jacohian of the treinsforHiatioiio This Justifies the name 
'’invariant" for the determinants h^^ and LeEour calls these the first 
generalization of the Darhoux invariants. 

Now for the equation (15) to admit an integral of the form 
z » u^, X "being an ar^bitrary function of x only, it is necessary 
and sufficient that the coefficients A of ( 20 ) are all zero. This is 
evident if we "but consider 



Since X is arbitrary, and this equation must "be identically zero, it 
follows that 



which proves the assertion. Thus we see that the vanishing of the invar- 
iants plays the same role for (15) as it does for the second order 
equation, 

LeRo\ix then considers the case when x is a multiple character- 
istic varia"ble, so that (15) has the form 




1x0 -2 . X- « 








where ♦ * ' designate the differential multipliers, e.g. 




The first su"bstitutioa we have already considered, namely 




ll6 



tiiat 



We attempt to deccmpose <^q(z) into differential factors, so 

Jy ®^r-2 Tj 



Then we define 

z ■* u ^ 2 

^ ° SF ^ ^ 

Zq »Of^U^ 3 i ^ ^ 

JF ^ Uq 1 ; 

z3-oC«<u ^ ~ ^,^-or ocp u 

« • • 

and iiaing these w@ are ahle to redxace the equation. If does' not 

have such a deca!J 5 )ositlon, we may still reduce the equations by a method 
which will he described briefly. The proof of the following is rather 
lengthy and detailed and does not bear repeating here. It may be found, 
however in LeRous's work previously noted (footnote 22). 

Let Uq# S’ linearly independent integrals 

of the equaticna (l8) which define the conditions on the coefficients for 
the existence of a solution of the form, of (l6) . We then define the 
functions 



1 u 2 
- o 



u 



u 

o 

’^l 



^ o 

IT 



ko= Ju 
^ u 3 












>0 d\ 
d\ 



W~ ~T2 



ay 






} • • e 



U7 



Set 



« ^i^l-2 




The eq.uatioa 

d 1 i 1 ... 3x3/u \ _. 

admits as solutio»s the fuaactloas «eoU^» Then the set of transfor- 

mations Is easily defined, "by setting 





• ^ 



z 



i 



u 

o ^1 



^ 1 / z \ 

Jy li Jy [% ) 




1 1 • e . 

° 1 ^2 



'i-T. 



5 I ^ 1 

Sj'kTi^y Ti:i 



• O • 






If there exists a particular integcul of the form of Exaer, the invar- 
iants 1 will heccme zero after a certain mmber of Iterations, and 
the chain of transformations will stopo 
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C. Following the work of Darhoux and LeBoux, Dinl made an extension 
■| 

of the Laplace method to the linear second order equation in ein arhltrazy 
number of Independent varlahleso Dinl considered the eq.uation 
n -V 2 



n 



(27) 






and. posed the question, can we transform (27), into the form 



( 28 ) 



n 



1=1 



ki ^ + M0 + 

^ X, 



n 



i=l 



oC, 






z +L* 4-1 = 0 



where 



U8 



n 



(29) 






a , 






^ z 4 "bz . 



aEd tlie fuaRstioEs ^ (2n + 2) functicsas 

to 1)© cliosea ar'bitrajpily? If we assuma that this has heea accaaplished. 



substitute (29) into (28), and equate coefficients with the left side of 

(27), w© obtain the (n4l)(n+2) equations 

2 



(30) ^i^'J ^ ^ i,j = l,...,n^ 

n 

(31) kil) + a^M 4^ kj ^ + o<'i ~ \ j i=l,...,n^ 

J-1 ^Xj ^ 



(32) 







L a 



K 

J 



for the (2n+2) unknowns. We obserre, however, that in (30), the k^, 
and Sj appear only as products ^i^j# ^'3. not all of them are zero if 
(27) is to be effectively of second order. Thus we may select > arbitrar- 
ily one k, say k^, and set k^ »1, without affecting these relations; 
and hence there are actually only 2n-l unknowns a^^ and k^, . If we consider 
first the equations (30), we not© that for ns2, these are 3 ©quations 
in 3 unknowns, and this system is determined. For n » 3; however, the 

system is overdetermined, and on© would expect to have (n-l)(n-2) rel- 

2 

at ions among the coefficients A^j so that only 2n-l of the equations (30) 
will be independent. 

W© see then, that befor’© we may even begin to seek apply a 

"Laplace method” to equation (27), the coefficients must satisfy 

(n.-l) ( n«^) conditions. Dlnl shorred that In general these conditions may 
2 

be stated as 





I 



r 
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(33) W - Wrs " ■/V - - Was , » ' 2.3, 

J 

where -1^ * il; stud these si^is are determiiied hy the siga of 

the radicals of the rocts of the quadratic equation. 



0 0 9 ^ 



Here Aj^ is considered to he the first noa-iraiiishljng coefficient 

t 

i=sl, oo.,n. If all the are zero, the equations (33) ar© 

satisfied identically for all i = l, 

The classification of the vai‘ious allowah.le types of equation 
(27) hecomes the next important consideration. We list here without 
proof sceae of the more impccrtant properties of (27) derived hy Dini, 
First, if n>2, and all coefficients of (27) are real valued, then 
for any pair of variahles (2^,Xg) for which the partial differential 
equation is not parabolic, the equation moat always he of the same type. 
That is, ignoriag all pairs of variahles for wdilch the equation is para- 
bolic, and considering all remaining pairs of variahles, the equation 
must always he elliptic, or always he hyperbolic type, with respect to>i^ch 
pair of variables. Another property, true even for complex-valued coef- 
ficient si- if 0, and the partial differential equation (27) is 

parabolic with respect to the pairs then it will 

also ho parabolic with respect to the pair (Xj,,x^) , These properties 
follow from the conditions (33) •• 

Dini then defines the following terminology. Equation (27) is 
said to be of parabolic type provided that it is parbolle with respect to 
every pair of variables (a:^,Xg)j if there exists at least one pair of 
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variables (x ) smh that (27) is not parabolie -wltli respect to ' 
(Xy,Xg) then (27) is said- to be of elliptic ( hyperbolic) type if it Is 
elliptic (hyperbolic) with respect to every pair of variables (x^^Xs) 
for which it is not of parabolic type; otherwise (27) is said to be of 
mixed typ@« Then Dini proves that if (27) is not of parabolic type^, 
then either there exist exactly two systems of values of k and a which 
satisfy (30);, or there exists aona;> while if (27) is of parabolic, type, 
there exist at most one system of values of k and a satisfying (30) . 

In the first case, on© obtains one, system of values from the other merely 
by Interchanging the k’s and the a's. Denoting these values by (k,a) 
and (a,k) respectively, w® may call them conjugate setso 

Mow our first thought in attempting to reduce (27) to the inte- 
gration of two partial differential equations of first order, is that 
perhaps (28) will reduce to a first order equation In 0 alone. In ad- 
dition to the (n-l) (n-2) conditions previously noted as necessary to 

2 

redw;e (27) to' (28) and (29), this requires n-l further candltians 
to insure that <^2.* ^2* ***’®^n* ^ These conditions arise fr^ 

the (n+l) equations (3I) and (32) idiich determine only the two un= 
kncrwns b and M, Assuming that the equation is not of parabolic type 
and hence/ without loss of generality, that it is not- parabolic with res- 
pect to the pair (x^^^gJ# Dini shows that the n-l eonditiaas can be fotrraulated 
as 

iSh) 

and 



M - Mb 



B— 0 



I 
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(35) 



All - A^ 



12 



A 



Is 



22 



2e 



J “• A 
2 2 



G » A 
8 s 



•“0 ^6 “3 f o o • ^ 



vb#r# A = 

s 






r-1 






aad B - 



n 

Z 

r-1 



k ^ 






W« i»y look upoa L and the C?C^ ^ (28) as playing the role 

of the inyariaats. Thus to say , that the coaditioas (3^) and (35) 
are satisfied, is to say that the Inyariants of (27) are all zero. 

*/ \ B. 

If we denote D (u) = M^ k o u ^ these inyariants laay he writtea 

J=1 ^ Jx 

as 



L ss N - D*(^)^ 

(36) 

°^i - G^- k^h - D*(a^) . 



In order to start the chain of oquations when these inyariants are not 
all zero, Dlni imposes conditions similar to those Imposed hy us on the 
third order equations of Section lY. For equation (27) to he cascaded, 
Dini requires that conditions (35) "be satisfied. That is to say, the 
inyariants are all zero, hut the inyarlant L is not. Then 

equation (28) may he solyed for z, in terms of © and its first 
partial deriyatiyes. This expression is suhstituted into (29) and a 
new second order equation, for © , is obtained. Not only will this 
eqtiation he of the earn© form as (27), hut the second order coefficients 
Aij are reproduced. Hence our new equation is 




i 
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n 



n 



(37) ^ A I© + 21 <3,' ^ + N'e + h'= 0 



^ X. 



The new coefficients are given by 



G = G + K - A-kL 
8 s s s 8 - a j 



(38) 



K 






K -l-Z\ - MLa -®) 



H * H(l) fE^ - 1 ,^)^ 



where K 



s 



V F 



\ ~"s # gi'»’en by 



H 




Since the relations (30) used to determine (k,a) for (27) involve 
only the coefficient s A.,, it imsediately follows that if there exists 
a set (k,a) which splits (27) into (28) and (29)^, the same set 
(k,a) will suffice to split (37) analogously. It will be necessary 

C 9 9 

however to coapute a new M h j> ^ 210 V set of invariants L ^ 

® c 

and • The conditions under which the ^ ^ aJ^© zero become 



11 


^21 


^1 


■ -^1 




12 


^2 


Gg 


- Ag 


•*0 0 ^ ® ” 3 ^ 0 ® 0 ^ n 0 


■X3 


■^20 


g’ 

s 


-Ae 





(39) 
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which are quite similar to the conditions (35) • If these conditions 

« » f • I 

are satisfied, and further L = W -Mb -B =0, then (37) reduces 
to two first order equations as desired. If L* is not ssro, however 
we must again require that conditions (39) be satisfied before we can 
iterate the process to obtain a third second order equation of the form 
(27). 

Dini’s very remarkable result, however, is the following. If the 

conditions (35) and the conditions (39) are all satisfied, then we 

need Impose no further condtions to iterate the process indefinitely. 

That is to say, if the invariants and ^ are all zero, then all 

(sl) 

invariants will be zero for any positive number m of iterations. 

This result becomes evident if we compute, say, G** by the rule given 

s 

9 tt tv 

in (38) for G^, substitute G^ - Gg " A^, and G^ - A^ for the 

corresponding terms in (39)# and employ the rule for evaluating a deter- 
minant when one of the columns is a sum, of two or more columns. Each of 
the resulting determinants in the s\im. must vanish as a direct consequence 
of conditions (35) and (39). 

In summary then, in order for the equation (27) to be cascaded 
as many times as necessary for the L invariant to vanish, there must be 
satisfied a total of 

(n-l)(n-2) -f 2 (n-2) — (n-|r3)(n-2 ) 

2 2 

conditions on the coefficients, The first (n-l) (n-2) conditions are 

2 

necessary to reduce (27) to two first oi'der equations, while the remain- 
ing 2 (n-2) are necessary to permit the cascading to continue. We observe 
that for n - 2, no conditions are necessary, and this is the original case 




% 



% 




i-a 




conaidered l>y Laplace, Legendire and Darboux 



(2l^) ^ 

D. Finally, Mention should "be uade of the work of Burgatti who 
considered the equation of elliptic type. 
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^ z 4. h ^ z 

^ a? 



+ cz= 0, 



Burgatti found that the two expressions 



(41) 



H = 
K s 








2c 



are inTariant relatiTe to the transfonnation z » X z*. Further, if 

H and K are hoth zero, equation (40) laay he reduced to 

2 * 2 » 

£ L. 4. O z - 0 , 

:h^2 2 ^ 

which is Laplace's equation; if H is zero hut K is not zero, then 
(40) takes the farm 

d^z* I z* c*z*g 0 . 

d x^ ^ ^ 

while if K is zero, hut H is not zero, (40) hecoraes 










= o 






where OC an*! ^ Bxe functions of x and y 
In the notation of Section V, let 
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then 



Dx * 2 



^ + a 



D ’ = 2 ^ + h 



and the Invariants (l^l) then ‘becctae 



H 

( 1 ^ 2 ) 



K 



Dj(a) - l4 (i) = (B^a) - c) - (Bx(l>) - ?) 

3. Z 2 

i ^ 

D^(a) ^ DyU) - 2C = (D^ia) - c)4 (Dy(h) - c) 

'2 ” 2 2 



Since the change of variahle and the vanishing of invariants does 
not reduce the eq.uation to two first order equations, no consideration 
can he given to cascading the equations of elliptic type utilizing the 
invariants E and K, The work of Dini shows us, however, that this 
equation can he cascaded if we adiait the possibility of ccmplex coeffi*» 
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